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Abstract 

We review the recent developments of virtual Compton scattering. We focus on the 
kinematical regimes which look the most promising. The threshold regime gives access to 
the generalized polarizabilities of the proton. The hard scattering regime allows tests of 
perturbative QCD predictions and of the valence quark wave function. The Bjorken regime 
is closely related to deep inelastic scattering. It gives access to the off-forward parton 
distributions which generalize the ordinary parton distributions and may shed a new light 
on the spin problem. For each regime we discuss the experimental perspectives and the role 
of the Bethe-Heitler background. 
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1 INTRODUCTION 



Virtual Compton Scattering generally refers to any process where two photons are involved 
and where at least one of them is virtual. In the following we shall use Virtual Compton 
Scattering (VCS) in a more restricted sense, that is the reaction where a space-like virtual 
photon is absorbed by a hadronic target which returns to its initial state by emitting one 
real photon. The reaction can be accessed through electroproduction of photons, (e,e'7). 
In practice the target will be the proton but one could extend our considerations to other 
targets provided the experiments, which are very difficult, become feasible. The experimen- 
tal difhculty has two origins. Firstly, the cross section is suppressed by a factor a ^ 1/137 
with respect to the purely elastic case. Secondly, the emission of a neutral pion which 
decays into two photons creates a physical background which may prevent the extraction 
of the VCS signal. When the target is a proton one can use kinematical conditions where 
its recoil momentum is large enough to be measured in a high resolution spectrometer, thus 
allowing to separate the VCS events from the Tr^events. Even in this favourable case the 
experiments are still delicate and have become possible only due to the advent of the new 
generation of electron accelerators. This is why, even though VCS has been put forward 
many times by theorists, the experiments are just beginning. 

VCS was first proposed as a possible test of Quantum Electrodynamics (QED) at small 
distance |Bjo58 but it is amusing to note that the first calculation | Ber58| of the VCS cross 
section was motivated by the experimental program of the Mark-Ill accelerator in Stanford, 
where the electroproduction of photons was considered as an embarassing background of 
the electroproduction of pions, the physically interesting reaction. Later the VCS process 



was recognized as interesting by itself and was studied by the same authors [Ber61| in the 
framework of dispersion relations. The interest for this reaction then disappeared for a 
while due to the lack of experimental data. Actually VCS is a small background! 

Aft er the discovery of approximate scaling in deep inelastic scattering, VCS was pro- 
posed [ Geo71 | to get further insight in the light cone singularity of cur rent co mmutators 
and later the process was analysed using the operator product expansion |Wat8S], providing 
tests of Quantum Chromodyna mics (Q CD) and of the models of hadrons. In the framework 
of the hard scattering picture Bro80| , the first, and till now unique, calculation Far90a | 
of the VCS amplitudes in the framework of Perturbative QCD (P QCD) was stimulated by 
physicists who where planning the first experiment of VCS | vBi89(| using the high luminosity 
apparatus of the (abandonncd) PEGASYS project. 

The VCS program came to live thanks to the advent of the CEBAF accelerator. The 
high luminosity and high resolution of the latter allowed to propose a convincing program 
of VCS [Ber93|. Due to the relatively low energy, the best physics case here turned out to 
be the threshold region where the so called generalized polarizabilities |Gui95| of the proton 
could be measured. It was then soon realised that a part of this threshold physics was in 
fact immediately feasible |d'Ho95| at the MAMI accelerator. The data taking took place 
in 1996 and at the time of writing are still under analysis. In the meanwhile the first few 



events of VCS were observed [vdB95|. 

The last important development is the recent proposal |Ji97a |Rad96a| that VCS in 
the Bjorken regime could be used to extract new structure functions which reduce to the 
ordinary parton distributions in the forward limit. The first moments of these functions are 
related to the hadronic form factors and the second moment gives access to the total spin 
carried by the quarks in the nucleon. This is an unexpected development which illustrates 
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the interplay between the existence of a performant accelerator and the apparitions of new 
and promising ideas. 




Figure 1: (a) FVCS amplitude, (b,c) BH amplitude. 
In the (e, e'7) reaction on the proton 

e + p e + p + "f , 

the final photon can be emitted either by the proton, giving access to the VCS process, or 
by the electrons, which is referred to as the Bethe-Heitler (BH) process. If we call T'^'^ 
the amplitude of the (e, e'7) reaction, we have 

rpee'i _ rpFVCS ^ rpBH 

where T^^^^ and corresponds respectively to graphs (a) and (b,c) of Fig.(|l|). The 

meaning of the superscript FVCS is that the corresponding amplitude contains the leptonic 
current while VCS refers only to the subprocess 

p + 7* — > p' + 7 . 

The BH amplitude can be calculated exactly in QED, if one knows the elastic form factors of 
the proton. Therefore it contains no new information on the structure. Unfortunately light 
particles such as electrons radiate much more than the heavy proton. So the BH process 
generally dominates or interferes strongly with the VCS process, and this may jeopardize 
the interpretation of the (e, e'7) reaction, independently of the experimental difficulty. The 
only way out of this problem is either to find kinematical regions where the BH process is 
suppressed or to have a very good theoretical control over the interference between the BH 
and the FVCS amplitudes. 

Assuming that this problem is fixed, another difficulty of the theoretical analysis is the 
complexity of the VCS amplitude itself. For a spin S target, the number of possible helicity 
states is 2 X (25* + 1) x 3 x 2 since the virtual and real photons have respectively 3 and 2 
degrees of freedom. Parity invariance reduces this number by a factor 1/2 so that for the 
proton one needs 12 complex functions to specify completely the VCS amplitude. These 
functions depend, for instance, on the virtuality of the initial photon (Q^), the scattering 
angle (^7*7) and the energy in the proton-photon center of mass (Vs)- Extracting the 
whole set of functions would involve such a huge amount of very difficult experiments that 
it looks definitely out of reach. 
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In view of the two difficult points discussed above, the safe attitude has been to restrict 
the use of the (e,e'7) reaction to cases where the situation is simple enough to allow a 
fruitful interpretation. Until now 3 cases have been identified: 



The low energy region, below the pion production threshold, where the VCS ampli- 
tudes can be parametrised by 6 generalized polarizabilities which are real functions 
of Q2. 

The hard scattering region, where ^/s is large with respect to the typical hadronic 
scale, and 0^.^ is close to 90°. The virtuality needs not to be large. There one 
can test the prediction of PQCD in the framework of the Hard Scattering Picture 
(HSP) pro8C |. In particular the phase of the amplitude which is predicted to be 



large [ |Far90a[ |, |K7o91 | can be accessed through the interference between the BH and 



the VCS processes. 

• The Bjorken region, where s and are large, / s finite and 9^*^ close to zero. 
This is often referred to as Deeply Virtual Compton Scattering (DVCS). In this limit 
the amplitude depends only on 4 parton distributions (called Off Forward Parton 
Distributions, OFPD's). 

In Section 2 we give our notations and conventions. We present the theoretical frame- 
work and collect the formulae that are used in the rest of the paper. Sections 3 presents 
the threshold regime. We have tried to propose a qualitative introduction to the physics 
of polarizabilities by considering the real photon as an applied external field. The low 
energy parametrization of the amplitude in terms of six generalized polarizabilities is then 
discussed in details and we propose an experimental program to determine them. In Sec- 
tion 4 we explain the interest of the hard scattering regime and we use the diquark model 
to delineate the experimental possibilities. We also discuss real Compton scattering be- 
cause the experimental prospects are better in this case. Section 5 deals with the Bjorken 
regime. The DVCS amplitude is computed in the handbag approximation for which the 
non perturbative input are the OFPD's. We use the same OFPD's to compute the meson 
production amplitude. Using a simple ansatz for the OFPD's we present some estimates of 
cross sections in kinematics corresponding to existing and future facilities. We conclude in 
Section 6 with a discussion of the open problems and perspectives. 
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Heaviside-Lorentz 


Gauss 


Electron mass 


rUg ~ 


nie ~ 


Proton mass 


m = .938GeV 


m = .938GeV 


Planck constant 


h=i 


h=i 


Light velocity 


c = 1 


c = 1 


Proton charge 


e 


ec 


QED fine structure constant 


aQED = e^/An = 1/137 


OQED = Gq = 1/137 



Table 1: Constants used in the paper. 





4-momentum 


Helicity 


Spin projection 


Incident electron 


k 


h 




Scatterred electron 


k' 


h' = h 




Target proton 


P 


hN 


a 


Recoil proton 


p' 


h'N 


a' 


Final photon 


q' 


X' 





Table 2: Notation for the 4-momenta and polarizations. 



2 PRELIMINARIES 

2.1 Units and constants 



The system of units used in this paper is the (natural) Heaviside-Lorentz one [ Jac75 |, 
but occasionally we shall also write some formulae in the Gauss system so as to fix some 
misunderstandings. Table (|l]) summarises our notations or values for the various constants 
of this paper. 



2.2 Kinematics 

Greek letters are used to note Lorentz indices. For the space components we use latin letters. 
Summation over repeated indices is implicit and the metric is 500 = 1, gu = —1- We have 
collected in Table (||) our notations for the 4-momenta and helicities or spin projections. 
We often neglect the mass of the electron with respect to its energy, which is permitted 
except when the angle between the electrons or between one electron and the final photon 
is exactly zero. These situations are irrelevant experimentally, so the approximation is safe. 
As a consequence of this approximation, the electron helicity is conserved, which simplifies 
a lot the formalism. For the proton the same would happen in the hard scattering regime 
at very large values of s, but in the threshold region the use of the helicity quantum number 
is not so compelling. In this case we shall use the (rest frame) spin projection because it 
leads to a simpler angular analysis. 

The virtual photon exchanged in the VCS process has 4-momentum q = k — k' , while it 
is q — q' in the BH process. In the following we shall note with a Roman letter the modulus 
of a 3-vector, e.g. q = |q| . A hat on a 3-vector will denote the direction of this vector. 

We shall often work in the C.M. frame defined by p + q = 0, and by default every non 
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Lorentz covariant quantity will refer to this particular frame. The unit vectors of this frame 
are noted [e (1) , e (2) , e (3)] and are such that the components of the 3-momenta are 

q = q(0, 0, 1) , 

q' ~ q' (sin0, 0, COS0) , 

k — k (sinacos(/), sinasini/), cosa) , 

k' ~ k' (sin a' cos 0, sin a' sin (/), cos a') , (1) 

and the scattering angle of the electron is 6*6 = a' — a. The invariant quantities which 
appear frequently in the following are: 

= -q-,t={ci-e()\s^{v^<lf . (2) 

and we have the relations : 

, s — rr? s — rr? — Q"^ s + ~[- Q"^ s + m? 

2.3 BH and FVCS amplitudes 

For definiteness we consider the case of an electron beam. For a positron beam one must 
change the relative sign of the BH and FVCS amplitudes. In the one photon exchange 
approximation and in Lorentz gauge, the BH amplitude is 

T^" = ^e';L^''u{p')T,{p',p)u{p) , (4) 

where the spinors are defined in App. The polarization of the final photon is e'{q',X') 
and the vertex F is defined as 

T-'iK^K) = F^{X)Y + iF2{X)a''P{K' - K)pl2m , 

X = (i^'-i^)2,Fi(0) = l, F2(0) = 1.79. (5) 

From QED one has the expression of the non symmetric lepton tensor 

L^^ = u{k') [r—r—i^ —7" + 7"^ ^7 7^ —7'^) u{k) , (6) 

\ j.(k' + q' ) — nie + te ^.{k — q' ) — rUe + i€ J 

which explicitely satisfies the gauge invariance conditions 

= L^- (g, _ ql) = 0. (7) 

The FVCS amplitude has a similar expression 



J.FVCS ^ S^e';H^'-'u{k')j,uik). (8) 



Note that in Eqs.(^, ||) we have extracted a factor for convenience. The hadronic tensor 
iJ'"' depends on the proton structure. At this point we do not want to make hypothesis 
which are specific of a model. We only assume that it is described by a strong interaction 
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Hamiltonian Hg which depends on some generic charged field and its momentum (t/JjTt). 
Under a gauge transformation they transform as 



ip exp(— iea(a;))'0, tt — > exp(ieQ;(x))7r, (9) 

with analogous transformation for {ip* ^ tt*). The distinction between Fermi and Bose fields 
is irrelevant here. For definiteness we assume canonical commutation rules 

['K{x,t),ij{x\t)]^ ~i5{x~x'). (10) 

The explicit form of Hs is not necessary for the moment. We simply assume the following 
form 

Hs = J dxHs (tt, V, d,^, tt*, d,r) , (H) 

and that Hs is invariant under the transformations (^) when a is constant. In general Tig 
does not depends on the gradient of the momenta but including such a dependence would 
not lead to any difficulty, apart the fact that the formulae would be longer. Also one could 
include non localities in Tig which are stronger than the first derivative. 

To introduce the coupling to the electromagnetic field one imposes that the Hamiltonian 
equations of motion 

^ = ilHs,^\ . (12) 

be invariant under the gauge transformations (|^). This is achieved by the minimal substi- 
tution 

diip — > diip + ieAi'>p, 
dii/j* dill)* -ieAiip*, 

Us 7^s + je(V'*7r* - V7r)Ao, (13) 

where the gauge transformation for the electromagnetic field is 

A^ A^ + d^a. (14) 

We need to compute T^^'-^'^only to lowest significant order in e. Therefore when we 
perform the substitution ( p^ ) in Eq.(ll) we only need to expand the Hamiltonian to order 
e^. This leads to the effective Hamiltonian 

/2 
dxej^.A'' + ^S^,A^'A'' , (15) 



where we have defined 



gij ^ ^.^.* ^^^S ^2 9'^'Hs ^ ^ ^ 

^mo ^ 5°'^ = 0, (16) 
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and h.c. denotes the hermitian conjugate. If we note the field strength F'^" — d^^A^ — d^A^^^ 
Eq.(]l5|) imphes the field equations 

d^.F"" = ef + e'S''''A^, . (17) 

Note that the contact (or seagull ) term 5^" has only space components. This is 
characteristic of the Hamiltonian formulation where the variables are the fields and their 
momenta. In the Lagrangian formalism the time components S^^ , S^^ of the contact term 
are not zero but they disappear when one builds the Hamiltonian. Here the use of the 
Hamiltonian formalism is compulsory because we must use perturbation theory without 
too specific hypotheses about the structure. So even if we had started from the Lagrangian 
we would have been forced to go to the Hamiltonian at some point. 

By computing the effect of V in Eq.(p^) in perturbation theory to lowest significant 
order in e one gets the following expression for H^^'^ 



i/- = /^xe-^^'-^<p1.-(x,0) ^,^^,_^^^^^ r(0) 

h':l,uu = <p'\s^''mp>- (18) 

This can also be written 

H^'' = -^j d'xe-^'^-^ <p'\T[fix),rm \p > , (19) 

with T the time ordering symbol, and since the equations of motion are by construction 
gauge invariant, H^''' satisfies the gauge invariance conditions 

In the CM. frame one has p = —q and p' = —q'. The hadronic tensor, which depends on 
(p, p') and the spin labels (cr, cr') or helicity labels (ft-, h!) will be noted H^^'^ {q' ,a' ^q,a) 
or H^'^{q' , h' , q, h) when necessary . 



2.4 The VCS amplitude 

It is convenient to expand the lepton current which appears in Eq.(^) on the basis of the 
polarization vector of the virtual photon which are defined in App. ^ Using the expressions 
for the helicity spinors (App. ^ one gets 



u{k',h)j''u{k,h) ^ n{h,X)e"{q,X), 

A=0,±1 



\ Q 



n{h,X) = ( -Ae-'^'^a(A/i) +(5(A,0)^\/2ej ^ 



a{z) = 



VI + e + 2zy/l - e 



V2 



(k + k') 



l\2 



(k + k'Y + q2 



l + 2-ijtan2 ej2 



(20) 
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where e is the usual hnear polarization parameter. It is equal to 1 when the virtual photon 
is linearly polarized in the lepton scattering plane. We define the VCS amplitude by 

T^CS^y^ ^) ^ £;*(A')i/^'^£.(A), (21) 

which allows us to write 

yWC5(y) n{h, A)T^C'5(y ^ (22) 

Thanks to the decoupling of the lepton current in Eq.(p2|) one can take the limit — > 
in T^C'-S'. In this limit the transverse part r^'^'^(A = ±1) coincides (up to a factor e^) with 
the amplitude for real Compton scattering. 



2.5 Tensor expansion of the VCS amplitude 

It is often useful to have a general expansion of the VCS amplitude defined by Eq.(^). A 
convenient, though not unique, tensor expansion in the CM. frame and in Lorentz gauge 
has been proposed in Gui95| . For the longitudinal part (A = 0) it has the form 



+ i[h[e'*.{q X g')o'-e(l) + 6^£"*.ga.e(2) + 6?je'*.((7 x q')a.e{'i)] , (23) 
and for the transverse one (A ~ ±1) 



f^-iTVCS ^ a*e"*.e + a"s'-.qe.q- 



It -'I* 



+ i[b\e'* .qe.{q x g') + b'le'* .{q x q')e.q']a.e{l) 
+ i[hle'*.e + b'^e"* .qe.q']ff.e{2) 

+ i[b\e'* .qe.{q xq')+ b'^e'* .{q X g ')£•<? V-e(3) , (24) 



where the normalisation factor M = 2-\/poPo is for future convenience and a are the Pauli 
matrices acting on the nucleon spin. The helicity labels A, A' have been omitted to simplify 
the equations. 

The 12 scattering coefficients (a', fo'i,...) are complex functions of (q, q', cos 6') which 
become real below the pion threshold due to time reversal invariance 

2.6 The ( e,e'7 ) cross sections 

In this section, we are giving the observables of the (e, e' 7 ) reaction. Due to the small 
cross sections the experiments will require such high luminosities that the use of polarized 
hydrogen targets seems excluded. Moreover detecting the polarization of the outgoing elec- 
tron or photon looks impossible as this can be done only through another electromagnetic 
process, thereby decreasing by a factor ~ 1/137 the already small counting rates. Therefore 
it looks reasonable to restrict our considerations to the case of a polarized beam and to the 
detection of the recoil proton polarization. 

Let (^^, ^^) the 4- vectors which describe the polarization state of the initial electron 
and final proton. By definition is orthogonal to the particle momentum (that is ^e-^ = 
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^p-p' = 0)) and, in the particle rest frame, ^ is twice the expectation value of the particle 
spin . We define the following density matrices 

f^.^An = ^u{p',a')^-^^^^u{p',^). (25) 
and the Lorentz invariant quantity 

aa' X'hh' 

in terms of which the invariant cross section is 

rf-r ,e) = wrS=^^w^'^'(^+ ^ - ^' - '^')-^(^^'^')- (27) 

32-\/(p.fc)^ - m^mi «o Po % 

Note that in the above expression, the states are supposed to be normalized as 

< p\p' >= {27rf2pod{p- p'). (28) 

After integration over the energy of the final photon and defining 

M = M{Q,£,P)+M{Q,-£,P) = 2M{Q,Q), (29) 

one gets for the unpolarized cross section 

da _ (27r)-^ 

dQ'^dsdtd(f)d(f)e ~ 128m2fc2^^^A(s, -Q2,m2) ' ^ ' 

where the azimuthal angle of the electron (f)^ is arbitrary and 

A{x, y,z) = \/x'^+y^ + z'^- 2xy - 2xz — 2yz. (31) 

The experiment is often performed by detecting the electron in coincidence with the 
proton and the (e, e' 7) event is tagged by a zero missing mass, {k+p — k' — p'Y = 0. The 
corresponding laboratory cross section is then 

da _ (27r)-5 k'j„6 s - rrf_j^ ^^^^ 



The next interesting observable is the single (electron) spin asymmetry (SSA) because 
it should be measurable with the existing highly polarized beams. We define it as 

^ Ai(e.o)-.M(-r,o) 
7W(e,o)+M(-e,o)' ^ ' 

In practice is proportional to k which corresponds to a longitudinal polarization of the 
electron. Then, because of time reversal invariance and up to corrections of order uqed, 
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one has A = below the pion production threshold (see discussion of Section Q). Also A 
vanishes by reflection symmetry if the electron and hadron planes coincide. 

Finally for a fixed electron polarisation one can try to measure the average polar- 
ization V of the recoiling proton . The components along of this double polarization 
observable are 



2.7 Virtual photon cross sections 

The (e, e'7) reaction is somewhat unique among the electroproduction reactions due to the 
presence of the BH process. However, there may be situations where this process can be 
neglected in which case it is interesting to present the cross sections in the usual manner 
that is in terms of virtual photon cross sections. To this aim we define the following helicity 
amphtudes |Kro96 



$1 = e^T^^'^i 1 1/2, 1 1/2) , 
$2 = e^T^'^^(-l - 1/2, 1 1/2), 
$3 = e^T^^^{-l 1/2, 1 1/2) , 
$4 = e^T^^^i 1 - 1/2, 1 1/2) , 



$5 = e2r^'^^( 1-1/2, 1-1/2), 
$e = e^T^'^^(-l 1/2, 1-1/2), 
$7 = e2r^"^^(-l - 1/2, 1-1/2), 
$8 = e^T^^^{ 1 1/2, 1 - 1/2) , 



$9^e2r^^^( 1 1/2, l/2)go/Q, 
$10 = e'T^^^(-l - 1/2, l/2)go/Q, 
$1, = e2j,^/C5(_i 1/2, 1/2) VQ, 

<I>i2 = e2j.^cS( i_i/2, l/2)qo/Q, (35) 

where T^'-^^ {\' , h'j^, A, /iat) is the amplitude defined in Eg. (pT|) . The factor qq/Q which 
multiplies the longitudinal amplitudes ($g to <i>i2) is a matter of convention and is consis- 
tent with the polarization vector of App.^. 
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The unpolarized (e, e'7) cross section where one retains only the VCS contribution 
reads then 

= 7;— Tt) — ; ^ £ — ; ^ £ cos 26 — ; h v 2e(l + e) cos — ; — , (36) 

L.dk',d6dt 2tt " \ dt dt ^ dt ^ ^ ' ^ dt ' 



where the virtual photon flux is 

_ aQED Kdab S - 



2^2 koAab2mQ^l-e) ' 
and the various virtual photon cross sections [| are defined as 



(37) 



^Note that only ar and cti, can be considered as cross sections. The other ones may be negative. 

i) The cross-section for transverse photons (which, at = 0, reduces to the unpolarized 
cross section da/dt for real Compton scattering) 

i=l 

ii) The cross-section for longitudinal photons 

4=9 

iii) The transverse-transverse interference term 

= -| sRe m $7 - $8 + $3 *5 - K *6] . (40) 

iv) The longitudinal-transverse interference term 
dcT T T c 

= Sfte m ($1 - $7) + <i>w ($2 + $8) + *n ($3 - $5) + $^2 (*4 + $6)] . (41) 

The two-body phase space factor c is given by 

lo7r(s — m^)^ 
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3 THE THRESHOLD REGIME AND THE POLARI- 
ZABILITIES 

3.1 Introduction 

The best way to understand the physics case of the VCS at threshold is to imagine that 
the final photon plays the role of an external applied field This may look strange 

but is possible provided the applied field has the adequate asymptotic conditions. Near the 
threshold, the energy of the final photon of VCS is small. So the applied field and its electric 
and magnetic fields are almost constant in time and space. So VCS at threshold can be 
interpreted as eleetron scattering by a target which is in constant electric and magnetic fields. 
The physics is exactly the same as if one were performing an elastic electron scattering 
experiment on a target placed between the plates of a capacitor or between the poles of a 
magnet. 

From elementary physics we all know what happens: under the influence of the applied 
field, the charge J° and the current density J inside the target get modified. For a weak 
field this modification is linear in the field and the coefficients of proportionality arc the so 
called polarizabilities. For a uniform isotropic medium only two are necessary: the electric 
and magnetic ones. In the general case one needs a tensor which moreover depends on the 
space point. If we note SJ'^{x) the modification of the current density we have, in the linear 
response approximation 

SJ^ix)^ J d'^yP'"'{x,y)Al'=\y). (43) 

The polarizability tensor Pf^" actually quantifies how the system adjusts its internal struc- 
ture to the applied field. The way to measure it is obviously to measure 5J'^{x) and we do 
this by electron scattering, exactly in the same way as we measure the equilibrium current 
J'^{x) of a free target. In the latter case we know that the experiment can be interpreted 
in term of the Fourier transform of J^{x), given by 

/"(g) = j (i4a;e«?-^jM(a;). (44) 

Similarly a VCS experiment measures the Fourier transform of dJ^{x), the current induced 

by the applied field. 

At this point we can see how the VCS generalizes the real Compton scattering. Though 
we always think of the zero energy limit, in practice this energy is of course never exactly 
zero. The corresponding non zero frequency of A"^^* induces a non uniform time dependence 
in 5J^{x) which then radiates an electromagnetic field. The latter, having a non uniform 
time dependence, decreases only as 1/r. So it can be detected at large distance. This 
is the final real photon of the Compton scattering. What we get is again the Fourier 
transform of 6J'^{x) but now computed along the real photon line, that is for = |g| ~ 
by energy conservation. This is to be contrasted with the case of VCS for which and q 
are independent. Moreover since the photon is transverse in the real Compton scattering, 
only the transverse component of 6J^(x) come into play. 

Now we come to the difficult point of the problem, at least for the experimentalists. In 
the above reasoning we have tacitely assumed that 5J'^{x) contains only the response of 
the internal degrees of freedom of the target. The reason is that it is the interesting part. 
However if the target has a global charge and/or a global magnetic moment, 5J'^{x) also 
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contains a trivial part due to the global motion under the influence of the applied field. If 
we put a proton in an electric field the first effect that we observe is that it moves as a 
whole. Similarly the magnetic field produces a precession of the magnetic moment. Since 
in a weak field the acceleration is not large, if we scatter an electron on this object we 
will learn no more than what we learn with a fixed target! This problem is absent when 
one studies the polarizability of a macroscopic sample because it can be fixed in space by 
appropriate means, which is not possible for the proton. This absence of restoring force 
explains why the trivial response due to the motion as a whole of the position and magnetic 
moment dominates over the response of the internal degrees of freedom. This is the physical 



origin of the low energy theorem Low58 for VCS. 

According to the above development, it is clear that to extract a physically interesting 
information from a VCS experiment one will have to subtract the trivial part which is due 
to the global motion. This part must be defined carefully and this will be done in a precise 
way through the low energy theorem. As stated above, this theorem originates from the 
fact that at low energy the response is dominated by the global motion of the proton in 
the external field. So one can ignore the internal structure to calculate this part of the 
response. All what we need are the parameters which control the motion, that is the mass, 
the charge, and the magnetic moment. Once the motion is known, we have to compute the 
amplitude for scattering an electron on this moving proton. For this we only need the elastic 
form factors and they are known. The calculation in this way is rather involved, mainly 
due to relativistic effects. This has been done only in the case of real Compton scattering 



|Gel54|. In the VCS case, one prefers to go directly to the quantum derivation where the 
use of gauge invariance and perturbation theory greatly simplifies the work. Nevertheless 
we think it is useful to keep in mind that the LET has a well identified classical origin. 

3.2 Electron scattering in an external field 

Let us put the above considerations on a more quantitative level. In the presence of the 
external field the total |^ current is, from Eq.(|l7|) 

J^^ej^ + e^S-^M^^*. (45) 

Note that we have set A = A®^* in the seagull term as the other part of the field would 
only give a contribution of order or more. 

To get the full current we need to evaluate the effect of A"^^^ on the motion. In quantum 
mechanics this amounts to compute the modification of the nucleon state by A"^^^ . The 
perturbation is given by Eq. (]l^) but clearly the seagull term would again produce an effect 
of order e^. So we can write the relevant Schroedinger equation in the form 



dt 

with 



i^\t>={Hs + V)\t>, (46) 



y = e y dxj^'Al^* = eV. (47) 
Note that we work in Schroedinger representation so that does not depend on time. 



that is the one that scatters the electron 
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To specify tlie boundary conditions we assume that the electron scattering happens at 
time t. The scattering takes place on a hadronic state |t, in > whose time evolution under 
the perturbation due to A^^* is such that 

\t, in >^ \Ni > wheni -oo. (48) 

Similarly the scattering leads to a state |t, out > such that 

\t, out |iV/ > wheni +00. (49) 

To make the boundary conditions consistent with the dynamics, we assume an adiabatic 
switching of the interaction by multiplying V by e*^* when computing \t, in >and by e~'* 
when computing |f, out > . 

Using standard time dependent perturbation theory we get to lowest order 

\t, in> = e-^^^* (l-t f di'e'^^*Ve^*'e-^^^*') |iV, > 



— 00 



|t, out > = e-'"'* dt'e'"^''Ve'*'e-'"''' j \Nf > . (50) 

At low energy we know that the target feels essentially a constant electric and magnetic 
fields [E, B). In this limiting situation the perturbation is time independent and, combining 
Eqs.(^, 1^) we get 

<t,out\J^{r}\t,in> = e''(^^-^-)*{<iV/|ej^(f) + e'^^"(r)^^"*(f) 

Ef - Hs + le Ei - Hs + le ) 

and we can identify the induced current as the part which goes like e^. 

In Eq. (|5j) we see that if we insert a complete set of intermediate states between V and 
the state corresponding to the nucleon itself is singular. This pathology is due to the 
fact that we have taken the strict limit of a constant field and is the quantum manifestation 
of the global motion discussed previously. Since this part contains no new information, we 
define the intrinsic induced current by excluding the nucleon from the intermediate states. 
Since the next possible state is a nucleon and a pion, the energy difference is then finite 
and we can set e = 0. Thus we get (at t = ) 



SJL.ir^ = e'{<Nf\S^'^{r)Ar{r)\Nf> 



< Nf\V\n >< n\j^'{r)\N,y < Nf\j^'{r)\n >< n\V\N,> 



— ^ Ef — En Ej — E',-, 



(52) 



To simplify a little we set Ei — Ef — m which amounts to neglect the recoil energy of 
the proton in the electron scattering event. This is only justified at small so in the full 
treatment later we shall not make this approximation. 

We now consider separately the effect of a constant electric or magnetic field. The gauge 
potential corresponding to a constant electric field is 

A',,t^-f.E, le.t = 0. (53) 
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Using Eq.(p7|) and remembering that the seagull has no time component we get 

I >< n\ 
En - m 



SJL.^,(r^ - E <^/l-^-^>;><-l^-^(-)l^-> + CO., (54) 



where we have defined the dipole moment operator 



d = J drrf{r). (55) 
For a constant magnetic field the gauge field is 



So we get 



A,.,t = -\?^B. (56) 



where the magnetic dipole operator is 

i^=\\ drfxj{r). (58) 

If we succeed in eliminating the trivial part of the response, then a low energy VCS ex- 
periment will allow to measure the Fourier transform of the induced currents of Eqs.(|5^, 

We can compute from Eqs.(^ ^ ) the induced dipole moments. After averaging over 
the nuclcon spin projection a we get 



•J ... , III' 

For the magnetic dipole we assume that the seagull term is as in a scalar field theory, that 
is S"-'(r) = —2S{i,j)j^{r), and we get 

Sfl = ]^ j drr y. J ^ {I3para+ Pdia)B, 

_ ^ I < jV,a|/i|n,a'> P 

2 p 



Note that in the above formulae the polarizabilities a and /3 are in Heaviside-Lorentz units. 
To get them in Gauss units, one must replace e by ec- 
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3.3 The low energy theorem 

Here we abandon the semi-quantitative picture and go back to the true reaction. As ex- 
plained above we need to separate the trivial part due to the global motion of the proton. 
However it is easy to see that this separation is not gauge invariant. By itself this is not a 
problem since only the full amplitude must be gauge invariant, but it is much more com- 
fortable from the theoretical point of view to deal with two parts which are independently 
gauge invariant. 

To this aim we start from the expression of the full VCS amplitude, Eq.(^8|), and instead 
of separating the on mass-shell nucleon contribution we write 

H = Hb + Hmb, (61) 
with Hg, the so called Born term, defined by 

i/^"- = uip')r^ip',p' + q') + ^ r-(p- + q',p)u{p) 

\P + h) ^ ^ 

+ u{p')V'{p\p~ q') ]-^ + "^ V^ip q\p)u{p). (62) 

(p — q') — 

Clearly Hb contains not only the nucleon contribution but also the pair excitation term. 
Moreover the vertex T is now evaluated for off mass shell values of one of its arguments 
because neither (p' + q')^ nor (p — q')^ are equal to m^. The advantage is that, if one uses 
the vertex decomposition of Eq. (|6^) , it is a simple excercise to show that 

i^r = ^r?- = 0, (63) 

which in turn implies the gauge invariance of the non Born term Hnb- 

It has been shown in ref.|Gui95| that H^b is a regular function of the 4- vector q'^, 



which amounts to say that H^b has a polynomial expansion of the form (this also holds 
for qf^ but we do not need it for the moment) 

iJj^^^ = a^^q) + b'^J{q)q'- + C,{q)q"^q'^ + ■■■ (64) 
Since the final photon is real we have q'^ = q' and therefore we can write 

i?j^''5 = a^'' + 6^:V" + 0(q''). (65) 
The gauge invariance of Hnb implies 

= q'^a^^ + 0(q'') - q'a"'' + g^a^ + 0(q''). (66) 

The coefficients a'^'^are constant with respect to q' . So averaging Eq.(|6^) over the direction 
q' on the one hand or multiplying Eq. (^6|) by q' and averaging again on the other hand, we 
get 

= q'a"" = q'a"' , (67) 
which implies that a^'' = 0. This low energy theorem (LET), which was first shown in 



rcf. |Low58 , tells us that the expansion of H^^b, the unknown part of the VCS amplitude, 
starts at order q' . Since, for finite the Born part starts at order 1/q' we conclude that 
the first two term of the VCS amplitude are known as soon as the form factors which define 
the vertex T are known. 
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3.4 Multipole expansion of Hj^b 



The LET defines in a precise way the separation between the trivial and non trivial part of 
he VCS amplitude. The non trivial part begins at order q' and is contained in H^qs- There 
is of course a contribution of order q' in Hb but it is exactly known and therefore can be 
subtracted, at least in principle. So what we need now is an adequate parametrisation of 
Hjsi B ■ For this we use the multipole expansion so as to take advantage of angular momentum 
and parity conservation. Such an expansion implies the choice of a frame and the natural 
one is the CM. frame. 

Because Hjsib is a regular function of q and g', wc know from elementary analysis that 

dq'H^^^^{q',q)Yi,{q')^{^f' , (68) 

where ~ means "at least of order" . A similar property holds for an average over q. Here 
1^/ stands for the spherical harmonic of rank (/', m') . Since we are interested only in the 
first power of q', we can restrict our considerations to the multipoles of rank I' = 1, which 
amounts to say that we are keeping only the dipole contribution of the outgoing photon. 

To factorize the angular and spin dependence of H'^''{q'a', qa) we introduce the reduced 
multipoles 

47rAA4^^'''^)^(q',q) 

y (-)'/'+"'+^+'' < i-d' ia|5s><L'M',L-A/|5s> 



aa'M'M 



J dq 'dq V; ip'L'M', q OiJ^B ('f'^'- 9'^)^. {pLM, q) , (69) 
where the basis vectors {pLM, q ) are de fined in Appendix M. The Clebsch-Gordan 



coefficients are the same as in Ref.[ Edm57 |. In Eq.(|69|), L (L ) represents the angular 
momentum of the initial (final) electromagnetic transition whereas S differentiates between 
the spin-flip (5* = 1) or non spin-flip {S — 0) transition at the nucleon side. 

The index p can take a priori 4 values: p = (charge), p = 1 (magnetic), p = 2 (electric), 
p = 3 (longitudinal), but gauge invariance relates the charge and longitudinal multipoles 
according to 

q,^(3L',pL)S^^,^(0L',pL)5 ^ 

So if we define the new basis 

W^^pLMA) = V^{pLMA), P=l,2, 

W^^(0,LAf,g) = y^(OLM,g) + — y^(3LM,g), (71) 

q 

which satisfy qfj,W'^''{pLM) = 0, p = 0, 1, 2, we have the following manifestly gauge invariant 
expansion of iJj^'^ 



AtiM 9p'p'W^{p'L'M',q')gppW''*{pLM,q) 

p'L'M'pLM 
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^(_)i/2+.'+L+M <L-^', Lo\Ss X L'M',L- M\Ss > ''^)^(q',q), (72) 

Ss 

where the sum over {p' , p) is now restricted to (0,1,2). The expansion (ff^ ) is vahd even if 
the two photons arc virtual but for the real final photon only the multipole p' = 1,2 can 
appear. This of course does not mean that the charge multipole does not exist: it is the 
contraction of the basis vector W with the transverse polarization vector of the photon e' 
which vanishes. 

The selections rules due to parity and angular momentum conservation are 

5* = 0,1, \L' - S\<L<L' + S, i-y+^' = , (73) 

and in addition one knowns that the electric and magnetic multipoles are zero for L = L' = 
0. 



3.5 Low energy behaviour and generalized polarizabilities 

Using Eq. (^) and the definition of the basis vectors it is clear that for small values of 
q' the multipoles behave as 

A similar relation holds for small values of q. This follows from the fact that in the charge 
and magnetic vectors the spherical harmonic has rank / — L. For the electric multipoles 
one has the same behaviour, but this does not follow from such a simple argument because 
in this case the spherical harmonic has rank I = L ±1. One has to use the Siegert relation 



Sie37| (see also Ref.[Gui95|) which is a consequence of gauge invariance and which states 



that 

^^^^^"'"^'(q',q) = -/^f i^^"f"'"^'(q',q) + 0(q'^'+^) , (75) 

with an analogous relation for the virtual photon. Since q"'^ = q' for the real photon, and 
(7° < q for the (space- like) virtual photon, we see that at small (q', q) all the multipoles 

behave at least as q'^ q^. As we are looking for the part of Hnb which goes like q' and 
using the selections rules we find that, a priori, this involves the following multipoles 

rr(ll,00)l rr(11.02)l „(11,22)1 r7-(ll,ll)0 r7-(ll,ll)l 

^NB ' ^NB ' ^NB ' ^NB ' ^NB ' 

„(21,01)0 r.(21,01)l r.(21,21)0 r7-(21,21)l r7-(21,12)l /-^n 

^NB ' ^NB ' ^NB ' ^NB ' ^NB ■ \'^) 

So the low energy (but arbitrary q) behaviour of the non Born VCS amplitude could be 
parametrized by 10 functions of q defined by 

Limit of -l-H^r^'P^'^^iq', q) when q' -> 0. (77) 
q' q^ 



with ( p'L',pL)S as in Eq.(Q. 

However this parametrization is not adequate. The first reason is that the (p' = 2, L' = 
1) multipoles actually behave as q' because of gauge invariance. It often happens that in 
model calculations gauge invariance is violated by some approximation which otherwise 
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would be innocuous. In this case the definition (^7|) may lead to a divergence. The second 
reason is more serious. In a VCS experiment the functions defined in Eq.(^7|) are measured 
by taking the limit q' ^ at fixed q. When one lets q — > one would like to recover the 
polarizabilities measured in real Compton scattering, but in the latter case the experiment 
realizes the limit (q' = q — *■ 0) and there is no guarantee that the two limits are the same. 
It is not difiicult to show | Gui95| that for the charge and magnetic multipoles the two limits 
are actually the same. This is again due to the fact that in those multipoles the rank of 
the spherical harmonic \s I — L or V — L' . By contrast, from Eq.([75|) one can understand 
that it will not be the case for the electric multipoles due to the factor (?°/q which goes 
to 1 in the Compton case w hile it obviously goes to in the VCS case. This problem 
was overcome in Ref. [Gui95| by the introduction of mixed multipoles which were neither 
electric nor longitudinal. However the findings of a recent work |Dre98 allow for a much 
more elegant solution. 

First one introduces the generalized polarizabilities (CP's) for the multipoles which are 
not problematic, that is 



p(ll,00)l(q) ^ 
p(ii.ii)5(q) ^ 



-H 



(11,00)1 
NB 



(q',q) 



p(ii,02)l(q) 



q'=0 



1 „(ll,02)l. , X 



q'=o 



1 „(1141)S 
q'q^^ 



(q',q) 



(78) 



q'=0 



Next one considers the multipoles where only the real photon is electric. In this case one 
can use Siegert relation to relate them to the corresponding charge multipoles which again 
are safe. So one defines 



p(01,01)S(q) ^ 



q'q 



(01,01)5 
NB 



(q',q) 



p(01,12)l(q) ^ 



'=0 



1 



(01,12) 



(q',q) 



, (79) 



q'=o 



which, using Eq. (|75|) , specifies the low energy behaviour of H^g^^'^^ and -ff^jj'^^'^. In 
Ref. [ Dre98| it has been shown that p(ii'00)i jg ^ot independent and can be eliminated 



according to 



p(ll,00)l(-q^ ^ ^^p(0l,0l)l 
90 



J_ 2p(ll,02)l 



(80) 



with (Jo the limit of when q' — s- 0, that is go 



The remaining problematic multipoles are Tf^jj'^^-*^, i/(2i,2i)o^ ^(2i,2i)i_ rpj^^ great 



value of the work of Rcfs.|Drc97|[Dre98| is that their authors have been able to express 
the low energy limit of these multipoles in term of the generalized polarizabilities already 
defined by Eqs.(^, |7|). Using their results, the definitions (Q^ and the relations (^,|^ 
one gets the following low energy parametrization of all the needed multipoles 



H 



(11,00)1 
NB 



(q',q) = VS^P^W^^)^^^2pinmi^^)\+oW'l 



(11,02)1^ 
'NB 

r(ll,ll)S^ 
'TVS 



i?r.'^^^"(q',q) 



H 



(21,01)5 
NB 



(q',q) 



1 

90 ■ ■ \/2' 

q,q2p(ll,02)l(q)^O(q,2^^ 

q'qP("'")^(q)+0(q''), 
-q'qV2P(°l^°l'^(q) + 0(q''), 



21 



rr(21,12)l. , x 


= -q'q' 




= q'qP( 




= -q'^o 




-4 



P("'")°(q)+0(q''), 

2q!p(ll,ll)l(q) _ V2q2p(0l>l2)l(qA ^ ^(^,2^^ ^^^^ 
Thus only 6 CP's, that is 

p(01,01)5(q)^ p(ll,ll)S(q)^ pdim)^^^^^ piOl^m^^^^^ (g2) 



are necessary to give the low energy behaviour of Hnb- Moreover it was shown in Ref . [ Dre98 1 
that one has the relations 

F("i^"i)1(0)= f("'")1(0) = 0. (83) 



The relations of Ref.| Dre98 | which allow to eliminate 4 of the GP's introduced in Ref. |]Gui95 | 



can be proved by combining the nucleon crossing symmetry with the charge conjugation 
invariance. These symmetries hold in any relativistic theory but to exploit them one needs 
a covariant parametrization of the amplitude. The reason is that nucleon crossing relates 
the reaction 

7*+iV(p)^7 + iV(p'), (84) 

to the reaction 

7* + IV(-p') 7 + N{-p) . (85) 

So to implement the symmetry one needs an expression where the 4-momentum formally 
appears while in the partial wave expansion only the 3 momentum appears. Though the 
correctness of the relations makes no doubt, their physical origin is not yet clear. In 
particular one sees on Eqs.(^l|) that a purely magnetic GP parametrizes a purely electric 
multipole. Since Lorentz transformations mix the electric and magnetic field, one may 
suspect that the origin of the relation is relativistic but a physical interpretation is still 
lacking. 

3.6 Generalized polarizabilities in models 

In a static model one neglects all the effects which go like the velocity of the initial or 
final nucleon. This allows to derive simple expressions for the GP's. Here we recall the 
expressions derived in Ref. ||Gui"95 1 for the most important ones 



p(o^^o,0(^) . ./Ij: 1 (86) 



[< iV|4(0)|X X X\d,iq)\N > + < N\dM\X >< X\d,{0)\N >] , 

[< 7V|/x,(0)|X X X\f^,{q)\N > + < N\f^,{q)\X X X|^,(0)|7V >] , 

(87) 
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where the generahzed dipole moments are defined as 

,3ji(qO -0 



d{q) 

M(q) 



dr- 



, dr 
2.1 qr 



/(f)f, 
qr 



where ji is the spherical Bessel function and we have omitted a possible seagull term to 
simplify. 

In the limit q ^ 0, if we compare with Eqs.(^ |6^) we find the relations with the usual 
polarizabilities 



p(01,01)0(-Q^ 



2 a 
3^ 



2aG 
34 ' 



(89) 



3 e 



G 



The generalized polarizabilities have been calculated in various more or less realistic 
nucleon structur e models . A first e stimate was proposed using the constituent quark model 
Gui95|| . In Rcfs. |Mct97a| |Met97b|] , calculations have been given in the linear sigma model. 



Calculations and a discussion of the polarizabilities in Chiral Perturbation Theory were 
given in Refs.|Hem97a| |Hem97b|. 



3.7 Low energy expansion of the scattering coefficients 

The calculation of the observables is easier and more transparent in terms of the tensor ex- 
pansion defined in Eqs.(23, 2^). The r elation between the scattering coefficients {a\ b\ . . .) 
and the multipoles has been derived [Gui95| using a generalization of the method used in 
pion electroproduction |Bcr67|. Combining the results of Rcf . | pui95 1 and the relations of 
Rcf . [ Dre98| one gets 



'■NB 



I _ qcost^-gp , 
w R 1 — z r o 



^^■1- q-goCOS0''^^'2 



hi 



'-NB 



^NB.3 



^NB,3 



f a ^" 
cos 

V q 


j "-NB 


^%B.l = ~ 


"TVS, 2 


"NB,1 


'^NBa 


q COS — qo 


ft'* 

'^NB,3 


q — go cos 9 



,P("l-"l)"(q)qq' + 0(q''), 



cos e - 1 ) P(01^01'l (q) qq' + 0(q'') , 



2 sin 9 \qo 

(f - 9/|p("'"^)^(q)) qq' + 0(q'^) , 



^ I cosy-^j P("'")"(q)qq' + 0(q'') 
3 f a!p(1141)l(q) _ V2q2p(01'12)l(q) \ ^ 0(^/2^ 



2 sin 9 \qo 

0(q'') , 
3 



. l-^cos0)p("^")l(q)qq' + O(q''), (90) 

2 sm 6' V 90 
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where NB indicates that we consider only the non Born part of the scattering coefficients 
and where Eg = Q^/('?oq)- 



3.8 Observables 

We now investigate how one can analyse the (e, e'7) observables in order to extract the 
GP's. We recall that the observables can be expressed in term of the Lorentz invariant 
quantity A4{^e,£,p) defined in Section 2. We consider it in the CM. frame and we choose as 
independent variables (q', q, e, 9, (f>). The tilde will denote the value of a dependent variable 
in the limit q' = 0. 

For the polarization observables we consider an electron beam in a pure helicity state 
h — ±1/2, which amounts to take = 2/ifc/r7ie up to terms of order that we can neglect, 
(this yields from Eq.@ = l/2meu{k,h) (1 + ■y'^2h)/2u{k,h) ). Below the pion 

threshold only the recoil polarization V, Eq.(^), is interesting. Since V.p' = one only 
needs the space part oiV and we consider its components along the basis [e (1) , e (2) , e (3)] 
defined in Section 2. We shall denote 



AM{h,i) 



M 



= 2hk/me, = e(i) 



M 



= 2ftfc/me, = -e(j) . (91) 



We assume that, at fixed (q, e, 0, 0), the experiment is able to determine A4 and AA4 in 
the form 



_^exp 



— 9^ + 




q'2 


q' 




AM 







AM'o'^P + O(q') 



(92) 



Due to the low energy theorem, the threshold coefficients Ai-2, -M^i, AAi^2, AA4-i are 
known. If we define 



^BH+Born 
^^BH+Born 



M 



BH+Bo 



^BH+Born 



-.'2 



^BH+Born 



O(q'), 



^_^BH+Bor„ ^_^BH+Born 



q 



/2 



q' 



AMl^ 



BH+Bo 



O(q') 



(93) 



where the index BH+Born means that the quantity is computed by keeping only the BH 
and Born contributions in the (e, e'7) amplitude, then we have the following constraints 
due to the LET : 



_^cxp ^ ^BH+Born^ AM!"P = AX^^+B"™, k^l,2 



(94) 



The information on the GPs is contained in A4q^^, and AAi'^^^ . These coefficients contain 
a part which comes from the (BH+Born) part of the amplitude and another one which is a 
linear combination of the GPs with coefficients determined by the kinematics. After some 
algebra Gui95 |,[ Vdh97a | one finds that the observables depend on the 6 following structure 
functions 
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|!^G^p(iMi)0(q) + |^GMP(Ol,Ol)l(q)^ 



^LT(q) 



(q) , 



(95) 



where Ge and Ga/ stand for Ge{Q^) and Gm{Q^). Clearly measuring the 6 structure 
functions defined in Eq.(p5|) amounts to measure the 6 independent CP's. For convenience 
we also introduce the combinations 



^LT — 



— 



G 



M 



TT 



2RGe 



TT 



RGi 
^Gm 

Gm 
RGe 

Gm ( ,QQp 



LL, 



{PIt - Plt) = -|Gm (^3P("'")^ + ^fpf"'")"^ , 

IGm ( 3^P("^")1 + ,/^«Ip("^")o^ 
2 \ qo V 2 q J 



(96) 



where R = 2m/Q. The observables ^ can t hen be written as (using the notation {x,y,z) 
for the directions (1,2,3) ) [|Gui95| , iVdh97a | 



■Mr --Mo 



BH+Born 



= 2K2 |fi [ePLL(q) - PTT(q)] + i^V2 - ^v^j ^2e {1 + e)PLT{q) 
AMl^'^ih, z) - AX|;'«+^°™(/i, z) 

= A{2h)K2[-vi^l^PTT{ci) + V2^j2e{\-e)PlT{q) + v^^2e [l - e)P'j;T{ci)] , 



The unpolarized cross section and the double polarizatio n observ ables yield 6 indepe ndent 



structure functions and not 7 as mentionned originally in Ref. [ Vdh97a | . Afterwards [Dre98|, four 



relations were found between the polarizabilities which yield then 6 independent observables to 
determine the 6 independent polarizabi lities. 

''Note the misprint in Eq.(4) of Ref. [ Vdh97a]: he should be replaced by e. 
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= A{2h)K2{vf^2e{l-e)PtT{<l) + v^^^l - P^^iq) 

+ «3Vl-£2^TT(q) + <v/2e(l-e)P;i(q)}, 

= A{2h)K2 {v\^2e (1 - e) Pij,(q) + ^f^l - ^TT(q) 

+ «|yT^F;-^(q) + t;4V2e(l-e)P;i,(q)}, 



with 



q 2m 

02 



Eg + m 



The angular dependent functions (wi, W2,- • • ) are given by 
vi = sin sin — kTCj cos 9 cos 0^ , 



V2 



Lo sin cos (f) — kxuj cos ( 



V3 ^ — [uj sin 6* cos 6* cos (/) — kyt^ (l — sin^ 6' cos^ 



sin cos (p (^Lu sin 6* — ky w cos 6* cos 

— sinS + kTW cos 6* cos 0, 
— cos 9 (^u sin 9 — kyw cos 9 cos 

kyw sin sin^ 0, 



(97) 
(98) 



and finally ^ 



Uj' = sin 9 sin sin 9 — hxuj cos cos ( 

V2 = k^ti-" cos sin 0, 

u| = kyti-" sin 0, 

i;| = — kTi^ sin sin cos , 









\ ,( I Ml 




\p.q' 


^ fc.g'). 


q'=o 


\k'.q' k.q' ) 


q'=0 



cjq — cj'A/k'^ — k|,, ky = 



2(1 



(99) 



(100) 



■'In Rcf.[Sui95| there is a misprint in Eq.(119) for uj, which should be corrected by replacing 
p' .q' by p.q . 
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3.9 Results for (e, e'7) observables in the threshold regime 



In the previous sections, we have outhned the observables of the p{e,e'p)j reaction below 
threshold and have detailed how the nucleon structure effect can be parametrized in terms 
of six independent polarizabilities. In this section, results will be shown for the p{e, e'p)^ 
observables in the threshold region in kinematics where the experiments are performed at 
MAMI [ |d'Ho95| and at CEBAF |Ber93[ . 

p (e, e' p)y : MAMI kinematics 
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Figure 2: p(e, e'p)7 differential cross section in MAMI kinematics : BH 
(dashed-dotted lines). Born (dashed lines) and BH + Born contributions (full 
lines) in plane (0 = 0'') and out-of-plane (0 = 60''). 



The experimental strategy of VCS in the threshold region consists of two steps. One 
measures the VCS cross section at several values of the outgoing photon energy. At low 
outgoing photon energies, the measurement of the VCS observables provides a test of the 
LET as can be seen from Eqs.(p2|, p^ . Once the LET is verified, the relative effect of the 
polarizabilities can be extracted using Eg. (p7|) . 

The predictions for the Bethe-Heitler (BH) and Born cross sections are shown in Fig. (||) 
for an outgoing photon of low energy. The BH cross section has a characteristic angular 
shape and displays two " spikes" . These " spikes" occur when the direction of the outgoing 
photon coincides with either the initial or final electron directions. In these regions, the 
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cross section is completely dominated by the BH contributions. In order to measure the 
VCS contribution, one clearly has to detect the photon in the half-plane opposite to the 
electron directions where the BH contamination is the smallest. 

The first absolute measurement of the VCS cross section on the nucleon performed at 



MAMI [ i'Ho95| has shown that radiative corrections provide an important contribution 
to the p{e, e'p)7 reaction. All diagrams that contribute to the cross section to order af„j 
have been calculated in Ref. [ Vdh98b|. A preliminary comparison to the MAMI data at 
low outgoing photon energy shows that the VCS cross sect ion can be understood by the 
radiatively corrected Bethe-Heitler (BH) + Born processes [ i'Ho97 | which provides a test 
of the LET. 



p (e, e' p)Y : MAMI KINEMATICS 
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Figure 3: Relative nucleon structure effect in VCS in MAMI kinematics. 
(^FULL represents the p{e,e'p)j differential cross section in the full model 
described in the text and ctbh+born represents the BH + Born contribution. 
The ratio calculated with all non-Born VCS diagrams is shown by the solid 
curve, whereas the ratio calculated with only the A diagram contribution to 
the non-Born VCS amplitude is shown by the dashed curve. 



For the sake of providing estimates of the nucleon structu re effect accessible through 
VCS, the relativistic effective Lagrangian model of Ref.| Vdh96 | is used here as it takes into 
account the nucleon structure effects to all orders in the outgoing photon energy. In this 
model, the p{e, e'p)j reaction below pion production threshold is described in terms of BH, 
Born and non-Born diagrams. For the non-Born diagrams, a tree-level calculation was per- 
formed consisting of exchanges of the resonances in the s- and u-channel and the exchanges 
of the TT*^ and a in the i-channel. For the resonances, the most important contributions 
come from the A(1232) and Z)i3(1520) for which the radiative couplings were taken from 
pion photo/electro production | Vdh95| ]. The model was tested against available data for 
real Compton scattering below pion threshold and yields for the proton polarizabilities at 
the photon point : ao — 7.29 • 10^*/m'^ and /?g = 1-56 • lO^^/m'^, which should be com- 
pared to the values extracted form experiment |MGi95| : a^^ = 12.1 ± 0.8 ± 0.5 • IO^^/to"^ 
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and f3'(^'' = 2.1 =F 0.8 T 0.5 • IQ-'^fm^. The underestimation of ac is probably due to 
the simplicity of the present tree-level effective Lagrangian model which docs not include 
non-resonant ttN intermediate states. 

To see the magnitude of the non-Born contributions, wc plot in Figs.(|) and (|) the 
BH-I-B-I-NB differential cross section {a full) relative to the BH-I-B differential cross sec- 
tion [gbh+b) for MAMI kinematics and for CEBAF kinematics respectively. In Figs.Q 
a stands for the fivefold differential p{e,e'p)j cross section da/{dV,^)Labd.\k l^ab {dflp)cM- 
The in-plane p(e, e'p)"f differential cross section ratios are shown for a momentum of the 
outgoing real photo n of = 120 MeV, which is close to but below the pion threshold. The 
MAMI experiment [|d'Ho95| was performed for a virtual photon momentum q = 600 MeV 
and a value of the virtual photon polarization around e = 0.61. In the MAMI kinematics, 
the non-Born contribution is a 10-15 % effect to the p(e, e'p)j reaction at backward angles 
(where the photon is emitted opposite to the elec tron d irections). 

The kinematics of the CEBAF experiment Ber93|, corresponds to a larger value of 
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Figure 4: Same ratio and curve conventions as in Fig. 
kinematics. 



but for CEBAF 



e w 0.95 due to the higher beam energy. The correspondingly larger virtual photon flux 
factor which goes like 1/(1 - e) increases the VCS contribution relative to the BH cross 
section. Thus one expects a larger effect in CEBAF kinematics as is also seen from the 
model calculations on Fig.(^) which are performed at q sa 1.06 GeV/c. 

If one wants to extract the polarizabilities from experiment, we have seen that an un- 
polarized experiment is not sufficient as it gives access to only 3 independent response 
functions. To separate the polarizabihties, one has to resort to double polarization ob- 
servables. Experimentally, at existing high duty cycle electron facilities with a polarized 
electron beam such as at MAMI, MIT-Bates and CEBAF, double polarization VCS exper- 
iments can be performed by measuring the recoil polarization of the outgoing nucleon with 
a focal plane polarimeter. 
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In Fig. (||) the double polarization asymmetry of Eq. ( p^ ) , where the recoil proton po- 



q = 600 MeV/c, q' = 120 MeV/c, e = 0.61 
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Figure 5: VCS double polarization asymmetry (polarized electron, recoil 
proton polarization along either the x- or ^-directions) in MAMI kinematics 
{Q^ = 0.33 GeV'^) as function of the CM angle between real and virtual 
photon. The BH + Born contribution is shown by the dashed lines. The result 
of the BH + Born + A contribution is shown by the dotted lines, whereas 
the BH + Born + D13 + a contribution is shown by the dashed-dotted lines. 
The total effect in the model is shown by the full lines. 



larization is measured either along the x- or z- directions, is shown at q = 600 MeV/c for 
in-plane kinematics. It is seen that the asymmetry where the final proton is polarized paral- 
lel to the virtual photon yields a large value (between 0.6 and 0.7). Remark that in Fig.(|^), 
the electric and magnetic contributions in the present effective Lagrangian calculation yield 
large and opposite effects to the asymmetry. 
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4 THE HARD SCATTERING REGIME 

4.1 Introduction 

Here we come to the hard part of this review. The threshold regime presented in Section 
H and the deep VCS which is the subject of Section [s] are going weU. In both cases the 
theoretical activity is important and while in the first case the experiments are running, 
in the second case they are already at the stage of experiment proposals. By an irony 
of history the threshold regime and the deep VCS should be considered as the wealthy 
descendants of the VCS in the hard scattering regime, a domain now somewhat sleeping. 
Though everybody agrees that this is a splendid physics case for VCS, it suffers from a severe 
drawback, namely the extreme difficulty of the experiments. It was the conclusion of the first 



feasibility study |Arv93| that no existing accelerator is adequate for such experiments. A 
new type of electron accelerator with a high energy (> 30 GeV), high duty factor 100%), 
very good energy resolution (much better than the pion mass) and high inten sity (> 100/^A) 



is actually necessary. Unfortunately no decision about the ELFE project | Arv95| , which 



meets with these requirements, is expected in the coming years. The consequence is that, 
despite its evident interest, the subject is in standby. In particular there is no motivation 
to develop the very complicated perturbative QCD calculations for this process. 

However we want to hope that one day the field will have the opportunity to take off and 
the role of this review may be to prepare the way. We shall do it by first introducing the 
physics of the hard scattering regime. Then we shall use the diquark model to get realistic 
estimates of the cross sections. We shall see that, even though the energy of CEBAF is 
too low to achieve a fruitful VCS program in the hard scattering regime, some preliminary 
studies can be performed with the polarized beam. Finally we shall switch to the case of 
real Compton scattering. The first reason is that the physics, though less rich, is closely 
related to the VCS because the relevant energy scales are (s, t, u) and not Q^, which 
can be zero. The second reason is that the experimental prospects are better in this case 



|d'Ho96||D'An97| 



4.2 Incoherence between the short and long-distance physics 

This regime is defined by requiring that all three variables (s, t, u) be large with respect 
to a typical hadronic scale, say 1 GeV. In this case there is a prejudice (actually proved 



in the case of elastic electron scattering |Stc97|) that the amplitude factorizes in a soft 



non-perturbative part, the distribution amplitude, and a hard scattering kernel which is 
calculable from perturbative QCD. Because of asymptotic freedom, the perturbative ap- 
proach must be to some degree relevant to the hard scattering regime. However, since the 
binding of the quarks and gluons in the hadrons is a long distance, non-perturbative effect, 
the description of the reaction requires a consistent analysis of both large and small scales. 
Here we follow the nice argument of Ref . pte97| . When the reaction is hard enough the rel- 
ative velocities of the participating particles are nearly lightlike. Time dilatation increases 
the lifetime of the quantum configurations which build the hadron. So the partonic content, 
as seen by the other particles, is frozen. Moreover, due to the apparent contraction of the 
hadron size, the time during which momentum can be exchanged is decreased. Therefore 
one expects a lack of coherence between the long-distance confining effects and the short 
distance reaction. This incoherence between the soft and hard physics is the origin of the 
factorization which is illustrated on Fig.(P). 
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Long distance Short distance Long distance 



Figure 6: Factorization of the scattering amphtude in the hard scattering 
regime. 

The calculation of the amplitude can be done using the Brodsky-Lepage formalism 
o80| which leads to the factorized expression 



T(A',/i'^,A, ft-Af) = / dxidyj (f>*N{yj) Th{X' ,h'^,yj, X,hN,Xi;s,t) 4>N{xi), (101) 



where (xi, yi) are respectively the momentum fractions of the quarks in the initial and 
final nucleon, Tu is the hard scattering kernel and 4>n is the distribution amplitude (D.A.). 
The evaluation of Eq.(lOl) requires a four-fold convolution integral since there are two 
constraint equations (xi +2:2 +X3 — 1 and yi +2/2 + J/3 = 1)- In Eq.(lOl) a sufficiently large 
momentum transfer is assumed so as to neglect the transverse momentum dependence in the 
hard scattering amplitude Tu- In this limit, the integration over the transverse momenta 
k^i (where k^i — 0) acts only on the valence wavefunction 

'^v{xi,X2,xz\ k^i, fc±2, fc±3) , (102) 
which is the amplitude of the three quark state in the Fock expansion of the proton: 



F >= *y I qqq > +'^gg \ qqq,qq > | qqq,g > 



(103) 



This valence wavefunction '^y integrated up to a scale fi (which se para tes the soft and hard 
parts of the wavefunction) defines the D.A. which appears in Eq.( |lOl| ) : 



(fk^i -^vixf^k^i) 



(104) 



For /i much larger than the average value of the transverse momentum in the proton, 
this function depends only weakly on fj, [ Bro80| and this dependence can be neglected. 

The interest of the formalism is that the distribution amplitude is universal, that is 
independent of the particular reaction considered. Several distribution amplitudes have 
been modeled using QCD sum rules | Che84 1 | ]Che89 1 [ Kin87[ . They have a characteristic 
shape and predict that in a proton, the u-quark with helicity along the proton hclicity 
carries about 2/3 of its longitudinal momentum (see Fig.(|^) ). 
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For the computation of the hard scattering amphtude Th, the leading order PQCD 
contribution corresponds to the exchange of the minimum number of gluons (in the present 
case two) between the three quarks as shown on Fig.(^. The number of diagrams grows 
rapidly with the number of elementary particles involved in the reaction (42 diagrams for 
the nucleon form factor, 336 diagrams in the case of real or virtual Compton scattering). 
Despite the large number of diagrams, the calculation of Th is a parameter free calculation 
once the scale Aqcd ~ 200 MeV in as{Q^) is given. Note that configurations with more 
than three valence quarks are a priori allowed but since this implies to exchange more hard 
gluons the corresponding contribution is suppressed by powers of 1/t. 




Figure 7: Model distribution amplitudes for the nucleon : KS (upper figure) 
and COZ (lower figure), as function of the valence quark momentum fractions 
xi and X2 {xi + X2 + xs = 1). 



There are two characteristic features of the Brodsky-Lepage model which are almost di- 



rect consequences of QCD: the dimensional counting rules [Bro73| and the conservation of 
hadronic helicities |Bro81 . The latter feature implies that any helicity flip amplitude is zero 
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and, hence, any single spin asymmetry too. The hehcity sum rule is a consequence of utiliz- 
ing the coUinear approximation and of dealing with (almost) massless quarks which conserve 
their helicities when interacting with gluons. Whereas the dimensional counting rules are 
in reasonable agreement with experiment, the helicity sum rule seems to be violated even 
at moderately large momentum transfers. The prevailing opinion is that these phenomena 
cannot be explained in terms of perturbative QCD (see, for example, Ref. [ 3iv89| ), rather 
they are generated by an interplay of perturbative and non- perturbative physics. 
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Figure 8: The hard scattering kernel in leading order PQCD. 



An interesting aspect of real and virtual Compton scattering is that they are the simplest 
processes in which the integrals over the longitudinal momentum fractions yield imaginary 
parts. The reason is that, as in any scattering process, there are kinematical regions 
where internal quarks and gluons can go on their mass shell. The appearance of imaginary 
parts to leading order in is a non-trivial prediction of PQCD which should be tested 
experimentally. As explained below, the (e, e'7) reaction with polarized incoming electrons 
seems to be a good candidate for this investigation. 



4.3 The diquark model 

The only calculation |Far90a | of VCS in the Brodsky-Lepage formalism is unfortunately 



not reliable due to a treatment of the singularities which leads to inaccurate results. This 
is discussed in details in Section 4.4. To get realistic predictions and to delineate what 
would be the most interesting VCS experiments we rely on the diquark model calculation of 
Ref. [ Kro96[ | . Due to the (relative) simplicity of this model it has been possible to analytically 



integrate the singularities, which eliminates the problems encountered in Ref. |Far90a|. 
In the diquark model Ans87| | baryons are viewed as made up of quarks and diquarks, the 



latter being treated as quasi-elementary constituents which partly survive medium hard col- 
lisions. The composite nature of the diquarks is taken into account by diquark form factors. 
Diquarks are an effective description of correlations in the wave functions and constitute a 
particular model for non-perturbative effects. The diquark model has been applied to a va- 
riety of processes and successfully confronted to data. Among these applications is a study 



of the nucleon's electromagnetic form factors | Jak93|. In fact, predictions are achieved for 
both the magnetic and the electric form factors. For the latter quantity no result is obtained 
in the pure quark hard scattering model because it requires helicity flips of the nucleon and. 
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Figure 9: Feynman graphs contributing to tlic VCS process j* p 7 p. 
Graplis with the two photons interchanged are not shown, and the thick line 
represents a diquark. 



in so far, the electric form factor also represents a polarization effect. In the diquark model 
helicity flips are generated through spin 1 diquarks. The diquark model is designed in such 
a way that it turns into the theoretically well established pure quark picture asymptotically. 
So the pure quark picture of Brodsky-Lepage and the diquark model do not oppose each 
other, they are not alternatives but rather complements. 

Representative Feynman graphs contributing to the gauge invariant hard-scattering am- 
plitudes are displayed in Fig.(||). The blobs appearing at the gD, jgD and 77!) vertices 
symbolize three-, four- and five-point functions. These n-point functions are evaluated for 
point-like diquarks and the results are multiplied with phenomenological vertex form factors 
which take into account the composite nature of the diquarks. The perturbative part of the 
model, i.e. the coupling of gluons (and photons) to diquarks follows standard prescriptions. 
The results for VCS shown below have been obtained with the same diquark parameters 
taken from previous applications of the diquark model |Kro96]. 

The first point to consider is the relative strength of the VCS and BH processes. A 
strong interference between them, though not uninteresting by itself, would obscure the 
interpretation of the results because it would prevent the familiar analysis in terms of the 
virtual photon cross sections defined in Eq.(|3^). For that purpose we plot in Fig.(p^) the 
difference between the full ep epj cross section and the VCS contribution to it divided 
by the full cross section. This represents the BH contamination. Clearly the dominance of 
the VCS contribution requires a high energy and requires to detect the outgoing photon in 
the half-plane opposite to the electron directions {(j) = 180°). 

The dependence of the cross sections is obviously the new information we can get 
from VCS by comparison with the real Compton scattering. In Fig. ([Ill) the dependence 
of the VCS cross sections is shown in a hard scattering situation (s = 10 GeV^, cos 9 = 
±0.6). The transition from real Compton scattering to large as predicted by the diquark 
model, is obviously non-trivial. In particular one observes that in the forward direction 
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Figure 10: The difference between tlie full photon electroproduction cross 
section and the VCS contribution to it over the full cross section vs. cosQ 
for Q"^ ~ \ GeV^, s = 10 GeV^, two values of the beam energy /cql, and for 
several values of the azimuthal angle 



(cosfl = 0.6), the variation is smooth, compatible with the diquark form factor. In the 
backward region (cos0 = -0.6) there is a rapid variation at small and the slope seems to 
be related to the scale of the diquark mass. So by studying the variation as a function 
of the angle, one can observe the interplay of the two characteristic scales of the model. 

The (L, LT, TT) cross sections have relative sizes with respect to the dominant trans- 
verse cross section which are not negligible and their dependence seems to merit a 
thorough investigation. However one must keep an eye on the very small value of the cross 
sections. Without a dedicated accelerator nobody will have the opportunity to contemplate 
the actual curves. Remark that at the real photon point, the diquark calculations predict 
a very small (non-zero) value for the photon asymmetry (TT/T) over most of the angular 
range which should be compared with the PQCD prediction shown in the next section. 

By contrast to the above comments where the BH is considered as a nuisance, the 
regions of strong BH contaminations offer an interesting possibility to measure the phases 
of the VCS amplitudes. These phases, due to the possible poles of the quarks, diquarks or 
gluon propagators, induce in the (e, e'7) amplitudes non-trivial phases beyond the phases 
due to the azimuthal angle dependence (see Eq. (^)). In other words, the perturbative 
phases manifest themselves in the fact that the T matrix is not self-adjoint. For the BH 
process on the other hand, T = obviously holds. As shown in |Kro96|, information 
on the absorptive part T — can be obtained from the electron single spin asymmetry. 
Since at CEBAF the electron beam is polarized a measurement of the asymmetry seems 
feasible. It is important to note that this possibility of measuring the absorptive part of 
the T matrix by the electron asymmetry follows from parity and time reversal invariance 
and does not rely on any model hypothesis. 

As is well known, a one-particle helicity state transforms under the combined parity and 
time reversal operation as 

|fc,A) — .r;(A)|A?,-A) , (105) 
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where 77(A) is ±1 depending on A, the spin of the particle and on its internal parity. So the 
combined parity and time reversal operations transform a given helicity state into a state 
with the same momentum but with reversed helicity. If the interaction is invariant under 
the parity and time reversal operations, the T-matrix elements for, say, a 2 — > 3 process 
satisfy the relation 

(fci, Ai; fcj, A2; fcg, AglTlfci, Ai; A;2, A2) 

(106) 



Let us now assume that particle 1 is a spin 1/2 one, say an electron and note cr(±) is 
the differential cross section for photon electroproduction with specified helicity of the 
incoming electron. Then, ignoring all kinematical variables, the cross section with particle 
1 in a definite helicity state is 

a(±)= J2 r'^'''(A'i,A^,A:,;±,A2)p. (107) 

The difference of these cross sections Act = ct(+) — ct(— ) may be written as 

ACT = 5Re Yl [T'''''^K■^'2,K■,+,^2)+l[v^T'■''^~K,~^2~K■,~,-^2r] 

{\i,X'-} i 

X r-^'^(A;, Ai, a:,; +, A2)* -Ylv^T^'^'n^X'^-X',, -Xi,; -, -A2)]. (108) 



If there would be no absorptive part of T, i. e., T — T\ then, according to Eq.(106), the 
difference of the two helicity cross sections and hence the electron asymmetry would be zero. 
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Figure 12: The electron asymmetry at CEBAF (6 GeV) as function of cj) for 
several values of cos 9 ior s — 5 GeV^ and Q"^ = 1 GeV^ . 

Therefore, the asymmetry measures the non-trivial phase as stated above. If we consider 
only the VCS contribution, we find (dropping irrelevant kinematical factors) 

Act ~ sin (f>^m ($i - $7) + $*o ($2 + ^s) + *n (^'3 - $5) + '^*i2 (*4 + ^e)] (109) 

Thus, we see that in this case Act measures the imaginary part of the longitudinal-transverse 
interference term whereas ctlt measures its real part. In other words. Act measures the 
relative phase between the longitudinal and transverse VCS helicity amplitudes. 

It turns out that the diquark model predicts very small values for Act in the pure VCS 
case. The ultimate reason is that the longitudinal helicity amplitudes are much smaller 
than the transverse ones for VCS as illustrated on Fig.([ll|). 

However Act can be strongly enhanced by the interference with the BH process. In the 
regions of strong BH contaminations it then essentially measures the relative phase between 
the transverse VCS amplitudes and the BH amplitudes. This enhancement can be seen on 
Fig.(|l|) where we show the electron asymmetry for a kinematical situation accessible in 
the future at CEBAF (s = 5 GeV^, /cql = 6 GeV, = 1 GeV^). As can be seen from 
that figure, the asymmetry is large in the diquark model for small values of | cos9\ (~ 0.2) 
and values of the azimuthal angle around 30°. The magnitude of the effect is of course 
sensitive to details of the model and, therefore, should not be taken literally. However, 
this suggests that preliminary tests of the theoretical predictions can be made with the low 
energy polarized beam of CEBAF even though the cross section is in this case dominated 
by the BH process. 

4.4 Compton scattering in PQCD 

The VCS estimates in the diquark model make clear that a measurement of VCS at high 
energy and large angle will be very difficult and will require a dedicated accelerator even 
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for moderate values of ( ~ 1 GeV^). At these rather low values of Q^, the cross section 
is dominated by its transverse part. Therefore, it is natural to consider real Compton scat- 
tering at high energy and large angle to study the hard scattering. From the experimental 
point of view, there are much better hopes to realize a real Compton scattering experiment. 
Firstly, the count rates are much higher than for electroproduction. Also, the required high 
intensity, high energy {E^ ~ 10 - 15 GeV) real photon beam might be feasible at the exist- 
ing HERA ring [d'Ho96||D'An97]. This motivated us |Vdh97b] to perform a leading order 
PQCD calculation for real Compton scattering. 

In a first stage, we simplified the calculations by approximating the x, y dependence in 
the gluon virtuality in as{Q^) by their average values for a given distribution amplitude. 
In this way we also avoid large contributions in the end-point region (x.^ w 0, w 1) where 
the asymptotic formula for as{Q^) is no longer valid. 

On ce th e hard scattering amplitude Th is evaluated, the four-fold convolution integral 
of Eq.(lOl) has to be performed to obtain the Compton helicity amplitudes. The numerical 
integration requires some care because the quark and/or gluon propagators can go on-shell 
which leads to (integrable) singularities. The different numerical implementations of these 
singularities are probably at the origin of the different results obtained in two previous 
calculations [ |Far90a[ |?ro9l| . In Refs. |Far90at|Far90b| , the propagator singularities were 
integrated by taking a finite value for the imaginary part -|-ie of the propagator. Then the 
behaviour of the result was studied by decreasing the value of e. To obtain convergence with 
a practical number of samples in the Monte Carlo integration performed in Refs.|Far90a| 
[ Far90b| , the smallest feasible value for e was e « 0.005. In Ref.|Kro91| the propagator 
singularities were integrated by decomposing the propagators into a principal value (off- 
shell) part and an on-shell part. To compare these methods, we implemented both of them 
and found for the +ie method differences of the order of 10% for every diagram as compared 
with the result of our final method. It is not surprising that, when summing hundreds of 
diagrams, an error of 10% on every diagram can easily be amplified due to the interference 
between the diagrams. To have confidence in the evaluation of the convolution of Eq.(lOl), 
we compared the principal value integration method with a third independent method. 
This third method starts from the observation that the diagrams can be classified into four 
categories depending upon the number of propagators which can develop singularities : in 
the present case this number is 0, 1, 2 or 3. Besides the trivial case of zero singularities 
which can be integrated immediately, the diagrams with one or two propagator singularities 
can be integrated by performing a contour integration in the complex plane for one of the 
four integrations. For the most difficult case of three propagator singularities, we found 
it possible to evaluate it by performing two contour integrations in the complex plane. In 
doing so, one achieves quite a fast convergence because the integrations along the real axis 
are replaced by integrations along semi-circles in the complex plane which are far from 
the propagator poles. We checked this method by also implementing the principal value 
integration method and found the same result up to 0.1% for each type of singularity. The 
principal value method was found to converge much slower and is more complicated to 
implement especially for the case with three singularities due to the fact that the three 
principal value integrals are coupled. 

Having exposed the PQCD calculational framework for Compton scattering, we now 
come to the calculations which are performed with several model distribution amplitudes 
denoted as CZ [ phe84[ , COZ |Che89| and KS [^inS^. 

The highest energy data which exist for real Compton scattering were taken around 
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Figure 13: Unpolarized Compton cross section on tlie proton for different 
nucleon D.A. : KS (full line), COZ (dashed-dotted line), CZ (dashed line) and 
asymptotic D.A. (dotted line). Data are from Refs.|Deu73||Dud83||Shu79|. 



5 GeV and are shown in Fig.dl^). Although the energy at which these experiments were 
performed is probably too low to justify a PQCD calculation, we nevertheless show the 
comparison with these existing data for illustrative purposes. 

Let us first mention that the hard scattering amplitude has the s-dependence (T ^ s~^) 
which leads to the QCD scaling laws [Bro73|, that is ^ ^ for Compton scattering or 
VCS. 

The unpolarized real Compton differential cross section (multiplied by the scaling factor 
s^) is shown as function of the photon CM. angle. We first remark that the result with the 
asymptotic D.A. (~ I2OX1CC2X3) is more than one decade below the results obtained with 
the QCD sum rules motivated amplitudes KS, COZ, CZ. The results with KS, COZ and CZ 
show a similar caracteristic angular dependence which is asymmetric around 90°. Note that 
in the forward and backward directions, which are dominated by diffractive mechanisms, 
a PQCD calculation is not reliable. Comparing the results obtained with KS, COZ and 
CZ, one notices that although these distribution amplitudes have nearly the same lowest 
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moments, they lead to differences of a factor of two in tlie Compton scattering cross section. 
Consequently, this observable is sensitive enough to distinguish between various distribution 
amplitudes. 
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Figure 14: Calculations for the polarized Compton cross section for two 
helicity states of the photon : A = +1 (thick lines), A' = — 1 (thin lines). 
Results are shown with KS (full lines), COZ (dashed-dotted lines), CZ (dashed 
lines). 



In Fig, (p^) , we show the polarized Compton cross sections for the two helicity states 
of the photon and for a target proton with positive helicity. We remark that for all three 
D.A. there is a marked difference both in magnitude and angular dependence between the 
cross sections for the two photon helicities. Consequently, the resulting photon asymmetry 
S defined as 



E 



^ (T,A = 1)-^(T,A = -1) 
(T,A = l) + f (T,A = -1)' 



dt 
da 
dt 



(110) 



where A is the helicity of the incoming photon and where t means a positive helicity for 
the hadron, changes sign for different values of 6cm as shown in Fig.(n3). This suggests 
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that the photon asymmetry might be a useful observable to distinguish between nucleon 
distribution amplitudes. 
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Figure 15: Calculations of the photon asymmetry for Compton scattering. 
Results are shown with KS (full lines), COZ (dashed-dotted Unes) and CZ 
(dashed lines) D.A.'s. 



The predicted sensitivity of the asymmetry to the nucleon D.A. can be use d in the 
extraction of a D.A. from Compton scattering data in the scaling region. In Ref.[ Vdh97b | 
we have outlined a procedure to extract a D.A. from Compton data in a model independent 
way by expanding the D.A. first in a set of basis functions and then using the angular 
informations of the cross sections to fit the expansion coefficients. It was seen that the 
precision on these coefficients is greatly improved when one measures both unpolarized 
cross sections and photon asymmetries. Such an experiment might be feasible at the HERA 
ring in the foreseeable future and might open up prospects to study the nucleon valence 
wave function in a direct way. 
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5 DEEPLY VIRTUAL COMPTON SCATTERING 



5.1 Introduction 

Deep-inelastic lepton scattering (DIS) in the Bjorken regim e (Q^ — > oo and xb — Q^/2p- q 
finite) is a powerful tool for the study of nucleon structure |3te95|. In particular, polarized 
DIS experiments have revealed that only about 23% of the nucleon spin is carried by the 
quarks |Sti96|. This has stimulated new investigations to understand the nucleon spin. 
Exclusive virtual Compton scattering in the Bjorken regime (DVCS) has been proposed 
recently | Ji97a||Ji97b||Rad96a| to access a new type of parton distributions, referred to as 
off-forward parton distributions (OFPD's) which are generalizations of the parton distribu- 
tions measured in DIS. Their first moments link them to the nucleon elastic form factors 



and it has been shown |Ji97a| that their second moment gives access to the contribution of 
the quark spin and quark orbital angular momentum to the nucleon spin. 

The off-forward or off-diagonal correlations of quark operators in the proton have a 
long history (see e.g. |Dit88||Bal88||Jai92] and references therein). They appear in the 
recent literature also under the name of nonforward parton distributions or double parton 



distributions |Rad97|. Despite their different names and definitions, the generalized parton 
distributions have in common that they include the additional degrees of freedom brought 
in by their non forward nature. 

It has been proposed that these OFPD's can also be accessed through the hard exclu- 
sive electroproduction of mesons (vr", p^,...) for which a QCD factorization proof was given 
recently [ Col97 |. This factorization is illustrated in Fig. (|l^b) and is valid for the leading 
power in Q a nd all logarithms. Furth ermore, it is valid f or all x b and thus generalizes pre- 
vious results I Rys93 1 1 Bro94 1 [ Rad96b | . According to Ref. |Col97], the factorization between 
hard and soft processes applies when the virtual photon is longitudinally polarized because 
in this case the end-point contributions in the meson wavefunction are power suppressed. 

In two very recent independent works |Col98||Ji98|, it has been shown that factorization 
also holds for the DVCS amplitude in QCD, up to power suppressed terms, to all orders in 
perturbation theory. It has furthermore been shown in these works that the factorization 
remains valid independent of the virtuality of the final photon, so that it also applies to the 
production of a real photon. 

In Ref. I Vdh98a | , the leading order of the 7, 7r° and Q leptoproduction amplitudes were 
calculated and first estimates for the cross sections were given using an educated guess for 
the OFPD's. These three reactions are highly complementary because apart from isospin 
factors, they depend on the same OFPD's. In the following, we present in more detail the 
7, 7r° and p^ amplitudes at the leading order in PQCD. We then explore the experimental 
opportunities which may allow to determine the OFPD's. 



^the index L means longitudinally polarized p° 
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5.2 Light cone dominance of the amplitude in the Bjorken Hmit 
and definition of oflF-forward parton distributions 





(c3) (c4) 

Figure 16: Direct (al) and crossed (a2) iiandbag diagrams for DVCS; (b) 
diagram for the factorized meson electroproduction amplitude; (c) leading 
order diagrams for the hard scattering part Th of the meson electroproduction 
amplitude. 



The leading twist contribution to the DVCS amplitude in the forward direction is given 



I Ji97b by the handbag diagrams shown in Figs.(|l6|)(al, a2). To calculate the corresponding 
amplitudes, it is convenient to use a frame where the virtual photon momentum and the 
average nucleon momentum (see Fig.(|l6|) for the kinematics) are coUinear and along the 
z-SLxis. Furthermore, in the Bjorken regime it is natural to express the momenta in terms 
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of the lightlike vectors 



= ^(1,0,0,1), 



(1,0,0,-1) 



where the light-cone components are defined by a± = l/^/2(aO ± 



momenta have the following decomposition : 



-2e r + ^ n'' 



where the variables m^, ^ and ^ are given by 



— 2 9 

m = m — 



A2 



2^ 



2e = 2e 




2Lb ' 
2 



(111) 

The physical 

(112) 
(113) 



(114) 



Remark that for a fast moving proton along the positive z-axis (such as e.g. is the case 
in an infinite momentum frame) the component will be very large whereas the P~ 
component will be negligible. For DVCS, the outgoing photon momentum is 



q^-A^'^-2(^ - r + ( ^ - C'Ti" ) - A': 



(115) 



In the Bjorken limit, the component of the virtual photon momentum is of order 
whereas the component 17+ is of order 1. Therefore, the operator product in the VCS 
tensor of Eq.([l9|) is dominated, as in DIS, by free quark currents separated by a lightlike 
distance. This is due on the one hand to the dominance of the 4-fold integral in Eq.(|Tg| ) by 
the region y+ ^ 1/1 9 I on the other hand to causality which forces > | E1177| . So, 
in the Bjorken regime Eq.(p^) reduces to a one dimensional integral along a lightlike line. 
DVCS and hard electroproduction of mesons in the Bjorken regime will therefore select the 
leading twist part of the matrix element of the bilocal quark operator represented by the 
lower blob in Fig, (p^) 



|- f dy-e^^p^y-{p\M~l)Ml)\p) 



y'^=y±=o 



41(7 



H'^{x,tt)u{ph+uip) + E'i[x,^,t)u{p)ia+^^u{p) 



2m 



+ (757 )a/3 



, A+ 

ff«(a:,C,i)u(p)7+75?/(p) -f Ei{x,^,t) u{p h^T^uip) 

2m 



, (116) 
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where ip is the quark field and u the nucleon spinor. In writing down the matrix element 
Eq.(116) we are implicitly working in the light-cone gauge n ■ A = ioi the gluon field A. 



If one works in a covariant gauge, one has to add a gauge link between the quark fields. 

The leading twist matrix element of Eq.( |116| ) is parametrized in terms of four OFPD's 
H'l , E"^ , H'^ , E'^ . These OFPD's are defined for each quark fiavor {q — u,d and s) and 
depend upon the variables a;,^ and t. The light-cone momentum fraction x is defined 
by k+ = xP+ and t = The support in x of the OFPD's is [—1, 1] and a negative 
momentum fraction corresponds to the antiquark contribution. Note that in writing down 
Eq. (|ll6| ), we have taken ^ = ^ to simplify the presentation. As shown in Eqs.(114), these 
two variables have the same value in the Bjorken limit. From Eq.(113) one realizes that 
^ is the longitudinal momentum transfer. As the longitudinal momentum fractions of the 
nucleons cannot be negative, ^ is bounded by 



< C < - — < 1 



m 



(117) 



A glance at Figs.(p^a, b) shows that the active quark with momentum k — A/2 has 
longitudinal {+ component) momentum fraction a; -l- ^, whereas the one with momentum 
k + A/2 has longitudinal momentmn fraction x — ^. As negative momentum fractions 
correspond to antiquarks, one can identify two regions according to whether |a;| > ^ or 
|a;| < ^. When x > ^, both quark propagators represent quarks, whereas for x < both 
represent antiquarks. In these regions, the OFPD's are the generalizations of the usual 
parton distributions from DIS. Actually, in the forward direction, the OFPD's H and H 
respectively reduce to the quark density distribution q{x) and quark helicity distribution 
Aq{x): 



H^x,^^Q,t = 0)^q{x) , 
H''ix,^^0,t = 0) = Aq{x) 



where q and Aq are defined as | Jaf96| 



y+=y±=o 



Aq{x) = |- / dy-e^-P^y'{pS^\\mi-ni54'{y)\pS\i) 



(118) 
(119) 



(120) 
(121) 



In Eqs.(12C, 121), p represents the initial nucleon momentum in the DIS process and S\\ is 
the longitudinal nucleon spin projection. From Eqs.(12C, 121) the quark distributions for 
negative momentum fractions are related to the antiquark distributions as q{—x) = — q(x) 
and Aq{—x) — +Aq{x). 

In the region — ^ < x < ^, one quark propagator represents a quark and the other one 
an antiquark. In this region, the OFPD's behave like a meson distribution amplitude. 
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5.3 Sum rules for off- forward parton distributions 



By integrating Eq.(116) over x one gets the following relations between the first moments 
of the OFPD's and the elastic form factors (for one quark flavor) 



J dxm{x,^,t) = F-'(i) , j dxE''{x,^,t) = F| 



(t) , (122) 



^ dxH'^ix,^,t) = g\{t) , r dxE'^{x,tt) - h\{t) , (123) 
-1 J-i 

where F^, are the Dirac and Pauli form factors respectively, g\ is the axial form f actor 
and is the induced pseudoscalar form factor. Note that in the sum rules of Eqs.(122) 
and (123), the ^-dependence drops out. 

According to the considered reaction, the proton OFPD's enter in different combinations 
due to the charges and isospin factors. For the DVCS on the proton, the combination is 

H^DVCsi^^U) - Ih-/p + ^H'^^P + ^H^/P , (124) 

and similarly for H, E and E. For the electroproduction of and -k^ on the proton, the 
isospin structure yields the combination 

ff^„(x,e,i) = ^-H-'P + \h'^'p'^ , Hl,{x,i,t) = i= i^-H^/P + \h'^/p 

(125) 

and similar for E and E. The elastic form factors for one quark flavor have to be related 
to the physical ones. Restricting oneself to the u, d and s quark flavors, this yields 

J^"/P = 2Ff + F^ + Fl , F^'P = 2F^ + Ff + F^ , (126) 

where Ff and i^" are the proton and neutron electromagnetic form factors respectively. 
The strange form factor is given by Fl^P = Ff , but since it is small and not so well known 
we set it to zero in the following. Relations similar to Eq.(126) hold also for the form factor 

2 ■ 



For the axial vector form factors one uses the isospin decomposition 



9l^' = + ^gl , gT = + ^gl ■ im 

The isovector axial form factor g^^ is known from experiment ((/^(O) ~ 1.26) and for the 
unknown isoscalar axial form factor g^ we us e the quark model relation : g^{t) = 3/5 gA{t)- 
For ft,^ we have relations similar to Eq.(127). 

The second moment of the OFPD's is relevant for the nucleon spin structure. It was 
shown in Ref . [ Ji97a | that there exists a (color) gauge- invariant decomposition of the nucleon 
spin: 

1 = Jg + Jg, (128) 

where Jg and Jg are respectively the total quark and gluon spin contributions to the nucleon 
spin. The second moment of the OFPD's gives 



1 

Jg = -J^ dxx[H''{x,^,t = 0)+E''{x,^,t = 0)] , (129) 
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and this relation is independent of ^. The total quark spin contribution Jq decomposes as 



Jo 



-AS + L, 



9 ' 



(130) 



where 1/2 AE and Lq are respectively the quark spin and quark orbital contributions to 
the nucleon spin. We recall that AE is measured through polarized DIS experiments. So if 
one can measure the OFPD's through exclusive hard electroproduction reactions, the sum 
rule of Eq. (|12S|) will determine the quark orbital contribution to the nucleon spin. 



5.4 DVCS amplitude 



In the following, we assume as in Refs.| Ji97a | | Ji97b | that a factorization between hard and 
soft processes is valid for DVCS. A formal factorization proof has been given very recently 
|Col98|| Ji98]. Using the parametrization of Eq.(116) for the bilocal quark operator, the 



leading order DVCS tensor H'^^'^q £,yijs (defined by Eq.(|19D) follows from the two handbag 
diagrams of Fig. dl^a) as 



^L.O. DVCS 



dx 



DVCS 



{x,£.,t) u(p)-f.nu{p) + E^jjy(jg{x,£„t) u{p')i(j 



+1 



dx 



^DVCS 



(x, t)N{p')^.nj5Nip) + E^ycsix, t)N{p)^^ 



2m 



A-r 

2m 



N{p) 



■Nip) 



,(131) 



with £oi23 = +1. Tests of the handbag approximation to the DVCS amplitude have been 
proposed in Rcf. [ Die97 ] and experiments are pr opos ed to check the hypothe sis [ Ber98 



The leading order DVCS amplitude of Eq.(131 



as given first in Refs.||Ji97a[ [|Ji97q| , 
is exactly gauge invariant with respect to the virtual photon, i.e. q^, H'^'^q jjyQg = 0. 
However gauge invariance is violated by the real photon except in the forward direction. 
In fact H'^'^Q jjyQg Aj^. This violation of gauge invariance is a higher twist effect of 
order compared to the leading order term H'^'^q j^y^g- So in the limit Q'^ oo it 

is innocuous but for actual experiments it matters. Actually for any cross section estimate 
one needs to choose a gauge and this explicit gauge dependence for 7^ is unpleasant. In 
the absence of a dynamical gauge invariant higher twist calculation of the DVCS amplitude, 
we propose to restore gauge invariance in a heuristic way based on physical considerations. 
We propose to write 



^DVCS 



^L.O.DVCS 



(a • q' 



' TjXlJ 

1x ^L.O.DVCS 



(132) 



where is a four-vector specified below. Obviously H'^y^g respects gauge invariance for 
both the virtual and the real photon : 



q^H'^^ 



DVCS 



Ifj. ^DVCS 



(133) 
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Furthermore, as 



)\ ^L.O. DVCS 



(134) 



the gauge restoring term gives zero in the forward direction (A_l = 0) which is natural. We 
choose a'^ = p'^ because p ■ q is of order Q^, which gives automatically a gauge restoring 
term of order 0(l/(5^). This choice for is furthermore motivated by the fact that 



in the derivation of the leading order am plitude of Eg . (131), only the components at 
the electromagnetic vertices are retained |tJi97a | | Ji97b |. The above arguments lead to the 
following fully gauge invariant DVCS amplitude: 



^DVCS 



rrfj.i> 

^L.O. DVCS 



ip-(l') 



L.O. DVCS 



(135) 



We will illustrate the influence of this gauge invariance prescription, by comparing DVCS 
observables calculated with Eq.(135) and with the leading order formula of Eq.(131). 



5.5 Hard meson electroproduction amplitudes 

For the electroproduction of vr'^ and p^, mesons ^ at large values of a;^ (i.e. in the valence 
region), the leading order amplitude is given by the diagrams of Fig. (16:). This amplitude 
was calculated in Ref.|Vdh98a| where the following gauge invariant expressions were found 
for the current operators 



^0 

PI 



dx 



dz 



01 

2x1 



HP {x,^,t)u{p )j.nu{p) + {x,^,t)u{p )«a" 

Pl Pl 



2m 



u{p) , (136) 



dz 



p^ 



2x\ 



+1 



dx 



H^o [x, t t)u{p )-f.nj5u{p) + E^o {x, t)u{p h5-^u{p) j , (137) 

where $p(z) and ^-niz) are the p and tt distribution amplitude (DA) respectively. For the 
pion, recent data | Gro98[ for the 77^7*7 transition form factor up to = 9 GcV^ support 
the asymptotic form: 

$,(z) = V2/,6z(l-z), (138) 

with /tt = 0.093 GeV from the pion weak decay. For the rho meson, the DA is not as well 
known. A recent theoretical analysis [ Bal96| updating the QCD sum rule analysis by 



^Note that although we will only show results for the and in this work, the formalism also 
applies for other vector mesons such as the uj and <jf> or other pseudoscalar mesons such as the 77. 
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including 0(q!s) radiative corrections, favors a DA that is rather close to its asymptotic 
form. In the calculations shown below, we use the asymptotic DA for the rho : 

$p(z) = ^/2/p6z(l - z) , (139) 

with fp = 0.153 GeV determined from the electromagnetic decay 



the CZ DA for the p |Che84 , which is more concentrated near the end points, would yield 
an amplitude larger by a factor 5/3. 

The ( e, e M ) cross section for a meson M is given by the formula of Eq.(^). The 
L,TT,TL cross sections appearing in Eq.(|3^) are calculated from the helicity ampli- 
tudes for meson electroproduction. In the following, we give only predictions for cr^ as, 
at large Q^, only the meson production by a longitudinal virtual photon survives and be- 
cause the factorization theorem [ ]Col97 only applies for a longitudinal virtual photon. The 



longitudinal j* + p —i- M + p two body cross section daL /dt of Eq. ( p6[ ) is 

^ = 1. , ^ J E E 1 1 ^^'(^ = 0) • = 0, 4; , (140) 

where ep(A = 0) is the polarization vector for a longitudinal photon (following the con- 



vention of App.A, Eq.(152)). It is expressed in terms of the lightlike vectors and 



|.^(A^O)^-i{..p^ + |^n^r (141) 



In Eq.(|l40|), J'^ is the current operator for production of a p° or 7r° defined in Eqs(136 



137). One can see from Eqs.(136, 137) that the large behaviour of the amplitude forp^ 



and tt" production is 



I e,{X = 0) • J^(Am - 0, h^- hM) - I e,{X = 0) • + j = --. (142) 

As s = + Q^(l — xb)/xb, this leads to a l/Q^ behavior of da^/dt at large Q^. 

5.6 Models for the off- forward parton distributions 

Ultimately one wants to extract the OFPD's from data but in order to evaluate elec- 
troproduction observables, we need a first guess for the OFPD's. We shall restrict our 
considerations to the near forward direction because this kinematical domain is the closest 
to inclusive DIS, which we want to use as a guide. In this domain, the contribution of E 
and E, which goes like A, is suppressed. So in the following we shall keep only H and H. 
Furthermore, we assume that, at small — t, H"^ is proportional to the unpolarized quark 
distribution for which we take the MRS (Sq) parametrization of Ref . | MRS93| 

xdv{x) = Ada;°-^^(1 - a;)4-5'^(l - 0.87V^-K 0.82a;) , 
X {uv{x) +dv{x)) = A,,dX°-^^{l-xf -^^{l + l4:A^+l6.99x) , 

xS{x) = 1.87(1 -a;)^"(l- 2. 21V^+ 6.22a;) . (143) 

The valence parts of the quark distributions are normalized as 



-1 r + l 

dxuv{x) = / dxdY{x) — 1 , (144) 
Jo 
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and the sea quark distributions are parametrized as us 
Our ansatz for the OFPD's H"^ thus yields 



2s, 



2s = 5/9 5. 



(145) 

and obviously the ansatz satisfies both Eq.(118) and the sum rule Eq.(122). 

Similarly, we assume that H'^(x,^,t) is proportional to the polarized quark distribu- 
tion. For the latter we neglect the not well known sea part and we adopt the valence 
parametrizations of Ref. |Gos97] 



Auv{x) = cos8d(x) 
cos 9£)(a;) 



uvix) ~ -^dv(x) 



Adv{x) — cos Qd{x) 



--dv{x) 



1 



(146) 



l + Hoil-x)^/ 

with Hq = 0.06 such that the Bjorken sum rule is satisfied. This yields the normalization 



/ dxAuvix) = 0.98 , / dxAdvix) ^ -0.27 , 
Jo Jo 

which leads to the correct proton to neutron magnetic moment ratio /ip//i„ 
ansatz for iJ' is thus 

H-/P{x, t) = Auv{x)gl/P{t)/gl/P{0) , H^/Pix, ^, t) = Adv{x)g''i^ {t) / g"/^ [0) , (148) 



(147) 
-1.47. Our 



and one can check that both Eq. (|ll9| ) and the sum rule o f Eq.( |l23|) are satisfied. Note 
that as suggested by the bag model estimate of Ref.| Ji97c |, we have omitted any depen- 
dence of the OFPD's on ^ which should be justified for moderate Q^. At high values 
of Q^, this approximation is expected to break down and the evolution has to be 
conside r ed. The evo l ution of the OFPD's has recently been studied by several authors 
I Rad97i iBaig^ [ |Bel98i |Fra98 1 . 
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5.7 Results for DVCS and hard meson leptoproduction observ- 
ables 

We first show results for the electroproduction cross sections. In Fig.(pT|), the total lon- 
gitudinal p1 electroproduction cross section is shown as function of the CM. energy W for 
different values of Q^. At high CM. energies, it is well known that the Perturbative Two 

Y* + p ^ Pl + p : = 6, 9, 17 GeV^ 




10 



W in^GeV 



Figure 17: Total longitudinal cross section for p° electroproduction. Data 
from NMC [|Arn94| (tria ngles) at = 5,5 (highest point), 8.8 and 16.9 (low- 
est poi nt) GeV2, E665 |Ada97i (black circles) at = 5.6 GeV^ and ZEUS 
Der95| (open circles) at — 8.8 (upper points) and 16.9 GeV^ (lower points). 
Calculations are shown at = 6 GeV^ (full hues), = 9 GeV^ (dashed 
hnes) and = 17 GeV^ (dashed-dotted lines). The curves which grow at 
high W correspond with gluon exchange whereas the curves which are peaked 
below W « 10 GeV correspond with quark exchange. The incoherent sum of 
both mechanisms is also shown. 

Gluon Exchange Mechanism (PTGEM) dominates as soon as > 6 GeV^ and this mech- 
anism implies a l/Q^ behavior of cr^ |Bro94 



'ra96 



with E g. (5) of Ref. |Fra96| is shown in Fig.(|T 



The PTGEM cross section calculated 
We used the CTEQ3L parametrization 
| Lai95| for the gluon distribution. The PTGEM explains well the fast increase at high 
energy of the cross section. Note that the data point shown on Fig. (p7|) are consistent with 
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the more recent ZEUS data [ |Cri97 |. 

However, the PTGEM substantially underestimates the data at lower energies (around 
W ^ 10 GeV), where the quark exchange mechanism (QEM) of Fig.(p^) is expected to 
contribute since a;^ is then in the valence region. In Fig.(pT[) we show our predictions using 
the OFPD's proposed in the previous section. The incoherent sum of both mechanisms 
is also indicated. Fig.(^^ provides a rather strong indication that the deviation from the 
PTGEM of the data at lower energies can be attributed to the onset of the QEM. The good 
agreement between the data in this region and our prediction gives us some confidence in 
our factorized ansatz for the OFPD's. 
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Figure 18: Comparison between BH (dotted lines), DVCS (dashed lines) and 
total 7 (full lines) in-plane cross sections for COMPASS kinematics. 



The QEM is studied further in Figs.(p^,(|T9|) where the angular dependence of the five- 
fold differential DVCS ^, tt" and p° muoproduction cross sections are shown for kinematics 
accessible at COMPASS. In our calculations we add coherently the BH and DVCS ampli- 
tudes. From a phenomenological point of view, it is clear that the best situation occurs 
when the BH process is negligible. For fixed and s, the only way to favor the DVCS 



The small differences between the results shown here and those of Ref. [ G!ui97 | are due to the 
use of different parton distributions. 
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over the BH is to increase the virtual photon flux and this amounts to increase the beam 
energy. According to our estimate, the unpolarized {1,1 7) cross section in the forward 
region is dominated by the BH process in the few GeV region. To get a clear dominance 
of the DVCS process one needs a beam energy in the 100 GeV range. This is illustrated 




Figure 19: Comparison between (full lines), 7r° (dashed lines) and total 7 
(dotted lines) muoproduction in-plane cross sections for COMPASS kinemat- 



in Fig.(p^ where we show separately the BH, the DVCS and the coherent cross section 
at ^ 200 GeV and ^ 5 GeV^. In the valence region {xb w 0.3), the DVCS cross 
section is more than one order of magnitude larger than the BH in the forward direction. 
The BH cross section increases if one goes to smaller xg and becomes comparable to the 
DVCS cross section around xg « 0.1. This is due to the fact that in the BH process the 
exchanged photon has 4-momentum (q — q') which gives a 1/i behaviour to the amplitude. 
The value of t in the forward direction {tmin) is very small for small xb values and it 
increases with xb- P^-t xb — 0.7, the value of —tmin in Fig. (|l8|) is around 1 GeV^, which 
explains the sharp drop of the BH amplitude compared to the DVCS one. The sharp rise 
of the BH process in the forward direction at small xb puts a limit on the region where the 
DVCS can be studied experimentally. At COMPASS kinematics, xb — OA seems to be the 
lower limit. Although the BH is not a limiting factor at high xs, one cannot go too close 
to = 1 in order to stay well above the resonance region. 

In Fig. (^9|) , the xb dependence of the p° j ""^ ^ind 7 fivefold differential muoproduction 



cross sections are compared in the same kinematics as in Fig. (18). The xb dependence 
of the , 7r° and DVCS cross sections is a direct reflection of the x b dependence of the 
OFPD's. In the valence region [xb ~ 0.3), the p^, cross section, which is sensitive to 
the unpolarized OFPD's, is about one order of magnitude larger than the tt" cross section 
which is sensitive to the polarized OFPD's. The DVCS cross section is sensitive to both 
unpolarized and polarized OFPD's but, due to the additional electromagnetic coupling, it 
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is also about one order of magnitude below the pj, cross section. Although the pj, cross 
section is the easiest to measure, the 7r° and DVCS cross sections seem large enough to 
encourage a study of the feasibility of the experiments. The three reactions are highly 
complementary due to their different dependence on the OFPD's. 

Before considering the extraction of the OFPD's from the data, it is compulsory to 
demonstrate that the scaling regime has been reached. In Fig.(pO[) we show the forward 
7* + p M + p cross sections as a function of Q^. This illustrates the lever arm one 
has to test the scaling behavior, as the maximum value of will be given by the count 
rate limit of an experiment. For the DVCS, the leading order amplitude is constant in 
Q and is predominantly transverse. Therefore, the leading order DVCS transverse cross 
section dax/dt shows a l/Q^ behavior. To test this scaling behavior, one needs of course a 
kinematical situation where the DVCS dominates over the BH. 

The leading order meson electroproduction amplitude was seen to behave as 1/Q and is 
predominantly longitudinal. Therefore, the leading order longitudinal cross section daL/dt 
for meson electroproduction shows a behavior. Note that eventhough the p° elec- 

troproduction is not as clean as the DVCS due to its dependence on the wavefunction, 
a test of the scaling behavior would nevertheless be meaningful. The reason is that 

an uncertainty in the meson wavefunction would influence mostly the normalization of the 
cross section, but not its behavior in Q. 

y' + p ^ M + p : (M = Pl, Tt", Y) 
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Although a high energy such as planned at COMPASS is preferable, one can try to 
undertake a preliminary study of the hard electroproduction reactions using the existing 
facilities such as HERMES or CEBAF, despite their low energy. To this end, we compare 
in Fig. (pi]), the p° i ""^ and 7 cross sections as function of the beam energy at a fixed = 
2 GeV^ and xb = 0.3. Going up in energy, the increasing virtual photon flux factor boosts 
the p° , 7r° leptoproduction cross sections and the DVCS part of the 7 leptoproduction cross 
section. For the 7 electroproduction cross section the BH process is hardly influenced by 
the beam energy and therefore overwhelms the DVCS cross section at low beam energies. 



1 + p ^ 1 + p + M (M = Pl, 7t", Y) 
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Figure 21: Comparison between (full lines), 7r° (dashed hues), DVCS 
(dotted hues), BH (thin dotted hues) and total 7 (dashed-dotted lines) lepto- 
production in-plane cross sections at = 2 GeV^, xb — 0.3 and for different 
beam energies : = 200 GeV (COMPASS), = 27 GeV (HERMES), 

E^- = 6 GeV (CEBAF). 

To make a preliminary exploration of the OFPD's in the few GeV range and to test the 
scaling, the measurement of the pi leptoproduction through its decay into charged pions 
seems the easiest from the experimental point of view as the count rates are the highest. 
In particular, in the — 1 - 5 GeV^ range, the pi cross section can be used to delineate 
the scaling region. 

For the 7 leptoproduction at low energies, we suggest that an exploration of the DVCS 
part might be possible if the beam is polarized. The electron single spin asymmetry (SSA) 
does not vanish out of plane and is only due to the interference between the BH process 
and the imaginary part of the DVCS amplitude. Therefore, even if the cross section is 
dominated by the BH process, the SSA is linear in the OFPD's. To illustrate the point, 
we show on Fig.(|2^) the unpolarized cross section for a 6 GeV beam and the SSA at an 
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azimuthal angle = 120*'. In the 0° - 5" region, a rather large asymmetry is predicted, 
even though the cross section is dominated by the BH process. Of course this does not 
exclude the risk that at this small Q^, higher twist effects might be important. 




Figure 22: DVCS at CEBAF for out-of-plane kinematics (</> = 120°). Upper 
part shows the differential cross section : DVCS (dashed line), BH (dotted 
line) and total 7 (full line). In the lower part, the electron single spin asym- 
metry is shown for the gaugeinvariant amplitude (full line) and for the lead- 
ing order non-gaugeinvariant amplitude calculated in radiative gauge (dashed 
line) and in Feynman gauge (dotted line). 



To illustrate the effect of the gauge restoring term in Eq.(135), we have plotted in 
Fig.(|2^) the SSA both for the gauge invariant and non gauge invariant amplitudes. For the 



non gauge invariant amplitude of Eq.(131), the SSA is shown both in the radiative gauge 
and in the Feynman gauge. As expected, all predictions are identical at small angle. At 
larger angles the gauge dependence clearly shows up, especially in the Feynman gauge. 

Finally, we recall that the real part of the BH-VCS interference can be accessed by 
reversing the charge of the lepton beam since this changes the relative sign the BH and 
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FVCS amplitudes. This leads to the following charge asymmetry 




Figure 23: e+e^ asymmetry at HERMES. Upper part shows the in-plane 
c+ differential cross section : BH + DVCS (fuU hne), DVCS (dashed line), 
BH (dotted line). Lower part : e+e^ charge asymmetry with same curve 
conventions as in upper part of figure. 
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6 OUTLOOK AND PERSPECTIVES 



As we have seen in the previous sections, virtual Compton scattering is a very versatile 
tool. 

At threshold it gives us access to a new kind of electron scattering experiments where the 
target is submitted to an external electric or magnetic field. The outcome of the experiments 
are the form factors of the currents induced by the applied field. This generalizes the notion 
of polarizability in the sense that the spatial distribution of the induced current can be 
measured while in real Compton scattering only the dipole moment is accessible. For the 
proton there are six CP's which are functions of Q^. Their experimental determination 
is now under way but this is a delicate task. As explained in Section ^, the response at 
threshold is dominated by the global motion of the proton in the external field. Scattering 
on this moving object teaches us nothing more than what we already know such as its mass, 
its charge and its elastic form factors. To extract the true information one has to subtract 
this trivial but dominating part and this supposes very accurate experiments. To obtain 
the CP's, with an accuracy of 15% one needs cross sections with an accuracy of 3%. 



The first experiment has been performed at MAMI [i'Ho9£]|Lhu97a| and at the time of 
writing, the analysis has not yet reached a stage where the polarizabilities can be extracted 
from the data. On the other hand the low energy theorem is well satisfied by the very low 
energy data |d'Ho97t . To our knowledge it is the first time that this theorem is tested. 
The next experiments will take place at CEBAF |Ber93| in 1998 where the higher beam 
energy is expected to help by allowing a better choice of the virtual photon polarization 
e. Finally an experiment at small has been accepted at Bates [Sha97] with the aim of 
testing Chiral perturbation theory predictions. 

At present no experiments with polarization have been proposed. Though polarized 
beams exist now at the high duty cycle electron accelerators, we have seen that one needs 
to measure the polarization of the recoiling proton. Due to the low efficiency of the po- 
larimeters, below 10%, one can foresee that this interesting physics will not start before 
the unpolarized program has achieved its task. We have demonstrated however that these 
recoil polarization observables provide three independent observables and complement the 
three unpolarized observables. Therefore, the measurement of these recoil polarization 
observables holds promise to disentangle the six independent CP's. 

The hard scattering regime of VCS discussed in Section ^ gives access to the most 
elementary part of the nuclcon structure, its 3 quark valence wave function. Thanks to 
the large energy-momentum transferred to the nucleon, the complex configurations which 
build the full confinement are eliminated from the amplitude. The price to pay is that the 



cross sections are extremely small and without a dedicated accelerator like ELFE |Arv93| 
the experiments look impossible. However a preliminary attempt will be performed at 



CEBAF |Bcr93|. The beam energy is definitively too low to access the hard scattering 
regime but the use of the polarized beam to select the BH-VCS interference will provide 
the first test of the diquark model prediction for the phase of the VCS amplitude. Finally 
we recall that the physics case of real Compton scattering is similar to the one of VCS. 
Since the experiments are much easier it may be a good strategy to try to develop first the 
real Compton scattering in the hard scattering regime. Such experiments might become 
possible with a high intensity, high energy {E^ ~ 15 GeV) photon beam at the existing 
HERA ring |d'Ho96[FA^ . 

VCS in the Bjorken regime or DVCS is the most recent development in the VCS field. 
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Because of its strong connection with the time honoured inclusive DIS, it has become 
rapidely popular. In this regime one can access the so called off-forward parton distributions 
which may shed a new light on the nucleon spin problem. We have used a crude factorized 
model of the OFPD's to perform the first evaluation of the cross sections including the BH 
amplitude. Gauge invariance with respect to the outgoing real photon has been enforced 
in a minimal way in the sense that the corrective term is explicitely a higher twist effect. 
Its effect with respect to a calculation in the radiative gauge is about 1%. 

To devise a convincing test of scaling of the DVCS amplitude and to access the OFPD's 
from this reaction one needs a very high energy in order that the BH process be negligible. 
The 100-200 GeV muon beam of COMPASS seems adequate, at least for not too small 
values of xb- 

However it is tempting to make use of the existing machines to open the way. The trick 
is to climb on the back of the BH amplitude using its interference with the DVCS amplitude. 
Even though the cross section may be completely dominated by the BH, the interference 
remains linear in the OFPD's that we are looking for. The first possibility is the single 
spin asymmetry which we have found to be large and which selects the imaginary part of 
the DVCS amplitude. A proposal of experiment is now in preparation [ Ber98t . The energy 
might be too low to reach the scaling regime but since the experiment does not imply a 
huge effort it is worth an attempt. The second possiblity is the beam charge asymmetry 
which is also large and selects the real part of the DVCS amplitude. In this case one does 
not need an out of plane experiment. This could be a good case for HERMES or SLAC 
and their higher energy as compared to CEBAF would be a serious advantage. 

The high energy forward meson leptoproduction is closely related to DVCS because 
the same OFPD's enter the amplitudes. Moreover one can select the OFPD's by selecting 
the meson. In this way, the p^, electroproduction gives access to the unpolarized OFPD's 
and the 7r° to the polarized OFPD's. One needs another input, the meson wave function, 
and therefore the information one gets from this process is less pure than in the DVCS 
case. However it is a useful complement and we have evaluated the cross sections using 
the (conservative) asymptotic wave function, which for the pion is able to describe the 
recent data for the 77^^7*7 transition form factor | Gro98 |. We found comfortable cross 
sections especially for leptoproduction. In the meson case the BH problem is absent 
and therefore the measurement of the electroproduction may provide a nice test to see 
if one can access the scaling region already at CEBAF. 

An important point which has not much been discussed in the text is the problem of 
radiative corrections, which is important in the case of experiments with electrons. All the 
Feynman diagrams contributing to the cross section of order o.'q^jj have been calculated 
|Vdh98b| using the dimensional regularisation both in the infrared and in the ultraviolet 
regions. The only approximation is that the radiative effects have been neglected for the 
proton. This approximation is justified by its large mass and has been tested to be numeri- 
cally negligible using a model [Lhu97t]. The difficulty of the calculation is that it cannot be 
done analytically. After renormalization and cancellation of the infrared divergences, there 
remain multidimensional integrals over the Feynman parameters that one must evaluate 
numerically. The hard point is that the functions to integrate have sometimes wild varia- 
tions, of the order i^beam/^^e in a range of the order rrie/ Ei^eam- It has been necessary to 
develop very accurate integration methods to take care of this [Lhu97b|. Another problem, 
which is particular to VCS, is that the integrals may also have poles corresponding to the 
on shell propagation of the intermediate photon. To cope with that it has been necessary 
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to extend the integration in the complex plane of the Feynman parameters. 

As pointed out in the introduction of this review, three interesting kinematical domains 
of VCS have been identified and discussed in some detail in this review. The threshold 
regime, the hard scattering regime and the DVCS have a well established physics case and 
are worth the experimental effort. This does not mean that VCS experiments cannot be 
done in other kinematical regimes. This is the case for the resonance region (a/s < 2 GeV 
which will also be explored at CEBAF in 1998 |Ber93|. In this region the advantage of 
VCS over other probes is not so clear even if one forgets about the BH problem. The point 
is that in the resonance region the photon channel is, by unitarity, strongly coupled to the 
meson production channels. So the information is dominated by meson production and 
rescattering. This may be interesting by itself but it is much easier to study it directly 
using meson electroproduction rather than with the very difficult (e, e'7) experiments. 

It seems now that, despite its experimental challenge, VCS has reached the status of 
a powerful and promising probe of nucleon structure. Nevertheless some work remains 
to be done to comfort this position. In particular, the higher twist effects need to be 
estimated and the model estimates of the OFPD's should be developped. Finally, since the 
OFFD's enter the amplitude through a convolution, one has to devise an efficient way to 
extract them from the data. We hope that this review will stimulate such work and that 
the experiments in the coming years will confirm our optimistic view of virtual Compton 
scattering. 
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A Polarization vectors and spinors 

The unit vectors e(l), e(2), e(3) of the CM. frame have been defined in the text. The 
polarization vectors of a photon with momentum k are 



e(fc, 1) = e(fc, 2) X fc, e{k, 2) - P^^, e{k, 3) = k. 

\k X e(l)| 



(150) 



and the hehcity states are 



£(/c,±l)=T^ e(fc,l) ± ie(fc,2) , e(fc, 0) = e(/c, 3). 
V2 L J 

The 4- vector polarization states in the Lorentz gauge are then defined as 



(151) 



£'^(±1) = 







k/fco 

m 



(152) 



The helicity spinors are 



u{k, h) = 



y/ko+m Xh{k) 
2/iVfco - rn Xh{k) 



(153) 



with 



Xl/2(fc) 



cos(a/2) 
e''"^ sin(a/2) 



X-i/2{k) 



-e-^*sm{a/2) 
cos(q;/2) 



(154) 



where (a, (f) are the polar and azimuthal angles of the direction k. The rest frame spin 
projection spinors are 



u{k, a) 



\/fco + m Xcr 



\/ko - m a ■ k Xa 



Xl/2 



X-l/2 



(155) 



These spinors are the positive energy solutions of the Dirac equation {■y.k — jn)u(k) = 
and the Dirac matrices are those of Ref. |Bjo64|. 
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B Vector basis 

The vector spherical harmonics are defined according to 

yimik) = Y^ilX, iM I LM)Yix{k)e{n) , with = 0, ±1 . (156) 

A, /J 

The magnetic, (transverse) electric and (electric) longitudinal vectors of the multipole ex- 
pansion are, respectively 

MLM{k) = y^Mik), 



CLM{k) = \j^f;^yLM{k)-^^^^yttl{k). (157) 
They have the following useful properties 

MLM{k) = , £LM{k) = -ik X MLM{k) , jCLM{k) = kYLM{k), (158) 

^/ L(L + 1) 

where is the angular momentum with respect to k. The 4-dimensional basis {pLM, k) 
is defined as 

V^{QLM,k) = {YLM{k),Q), 
F''(lLM,fc) = (0,A4LM(fc)), 
V''{2LM,k) = {0,£LM{k)), 

V^{2,LM,k) = (0,£LM(fc)). (159) 
and one has the closure relation 

3 

Y,Y.9pp'^''ipLM,k)V''*{pLM,k') = g^-6{k - k') , (160) 

and the orthogonality relation 

j dkVi'*{pLM, k)V^{p'L'M', k) = Qpp'dLL'SMM' . (161) 
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